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SUMMARY 


In this study we have achieved the following: 

(i) Theoretical Aspect 

Using the plasma theory, we study the interaction between the high 
frequency waves and magnetohydrodynamic (MHD) in which a set of coupling 
equations resulted. On the basis of this formalism, we examined the 
modulation instabilities by an electromagnetic soliton in a current sheet 
and showed that there is a resistive instability, which eventually turns 
into an eruptive instability at the onset of the magnetic field 
reconnection. This mechanism could be used to explain the onset of a solar 
flare. 

(ii) Experimental Aspect 

In order to have better measurements for vector magnetic fields with 
high resolution on the solar surface, we have examined the possibility of 
increasing the present state-of-the-art optics so we may use it to improve 
the design and fabrication of a new space-borne solar vector magnetograph 
as part of SAMEX (Solar Active Measurements Experiments) . The detailed 
study is presented in Section 3 of this report. 

INTRODUCTION 

The operations of satellites, space shuttle and space station for 
communications, tracking, and surveillance can be interrupted, degraded, 
or even endangered as a result of powerful explosions, that we know as solar 
flares, on the surface of the Sun. These spectacular eruptions release 
shock waves, hot plasma clouds, highly accelerated atomic nuclei, and 
bursts of x-rays, untraviolet and visible-band electromagnetic radiation 
into interplanetary space. When the path of propagation of these high- 
energy emissions intersects our terrestrial environment it is impacted in 
various ways that may produce deleterious effects on the operations of the 
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systems mentioned above. 

Observations have shown that physical conditions in the solar 
atmosphere are strongly controled by the solar magnetic fields. The 
appearance of solar flares, seen in enhanced emissions in H-alpha and 
different lines in the ultraviolet and extreme-ultra-violet as well as in 
white light observations, provides indications of prevalent nature and 
importance of solar magnetic fields. Consequently, to understand the 
physics of active regions, the storage and release of flare energy and 
formtion of hot plasmas and mass ejections, it is imperative that we 
understand and study the evolution of the Sun's magnetic field. To achieve 
such as goal , we took both theoretical and experimental approaches . In the 
theoretical aspect, we investigate the triggering mechanism of the solar 
flares by studying the resistive and eruptive instability. This result is 
presented in Section 2 of this report. In the experimental aspect, we 
investigate the possibility of improving the resolution of the optics for 
the space-borne solar magnetograph. The results are included in Section 
3. 

2. THEORETICAL STUDY: Resistive and Eruptive Instability by 
Pondermotive Force with High-Frequency Plasma Oscillation. 

We investigated the subtle interaction between the magnetohydro- 
dynamic (MHD) and high frequency plasma waves and we derived the coupling 
equations for these phenomena. On the basis of this formalism, we examined 
the modulation instabilities by an electromagnetic soliton in a current 
sheet and showed that there is a resistive instability for constant-/ 
approximation (i.e., uniform magnetic field configuration), which 
eventually turns into an eruptive instability at the onset of the magnetic 
field reconnection. 
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2 . 1 Introduction 

It is known that the electromagnetic radiation from a plasma often is 
trapped by the self-induced electromagnetic field and forms a soliton 
(Krall and Trivepiece, 1978) . It is also realized that the width of a two- 
dimensional soliton is inversely proportional to its own strength. 
Therefore, the existence of a strong electromagnetic soliton expected in 
the thin currrent sheets (Li, 1985) • Under these circumstances, the 
current and solitons through the ponderomotive force, within the current 
sheet will be intimately coupled; these interactions will cause the 
instability of the magnetic field configuration in the current sheet. 
Therefore, the electromagnetic soliton-induced instabilty within a 
current sheet becomes an interesting and fundamental problem in space and 
laboratory plasma. 

In Section 2.2 of this investigation, we present the fundamental 
coupling equations between magnetohydrodynamics (MHD) and high frequency 
plasma and also discuss the fast and slow-time scale. In Section 2.3, the 
electromagnetic soliton-induced resistive instability and in Section 2.4 
we present the soliton-induced eruptive instability. We present a 
numerical example concerning solar coronal plasma in Section 2 . 5 and 
concluding remarks in Section 2.6. 

2.2 Coupling Equations Between MHD and High Frequency Plasma 
.he following equations govern a two component magnetized plamsa in a 
gravitational field with two time-scale approximation. 

8n e 

+ 7 • (n e v e ) = 0 , 

at 

^ n e -» 

+ 7 • (njVi) = 0 , 

at 


( 2 . 1 ) 

( 2 . 2 ) 
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3v e . _ e _ 1 

+ (v e • v)v e = — (E + - v e x B) - 
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. 7P e 
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+ K e i ( v e - v i) » 
(2.3) 
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+ (Vj • 7) Vi = (E + - Vi X B) “ 
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7Pi 


m i n i 
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+ K ei (Vj - V e ) 

1 3B 

7 X E = t 

c at 

_ 1 3E 4tt 

7 x B = — — + — (®n e v e - eriiVi), 
c at c 

7 • B = 0 


(2.4) 

(2.5) 

( 2 . 6 ) 
(2.7) 


Note these symbols and their meanings: n = the number density, v = the 
velocity, P = the thermal pressure, E = the electric field, B = the magnetic 
field, g = the gravitational field, e = the electric charge, m = the mass of 
of particle and v e t = collisional frequency between electron and ion. The 
subscripts e and i indicate the electron and the ion respectively. 

Based on the two time-scale approximation (Li, 1985) , all the 
quantities could be defined by 


A = (n e , nj; v e , v^- p e . Pi,* E, B) 
= A f + A s 


( 2 . 8 ) 


It could be assumed that the assemble average value of fast time-scale 
components over the slow time-scale vanishes; 


<A f ^ 0 . 


(2.9) 

Hence, in the slow time-scale region, the charge neutrality cond- 
ition holds, 

<n,e - n,e> = 0 <2- 10 > 

Thus, 
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n 


e / 8 


= n i/8 E n, 


( 2 . 11 ) 

where, the subscript s indicates the slow time-scale. 

In the meantime, we shall employ this concept to analyze the set of 
governing equations (2.1) through (2.7). Beginning at the electron 
component of plasma, Eq. (2.1) the equation becomes, 


9 

— (n s + n f ) + v • 

at 


(n s + n f ) (v 8 + Vf ) 


= 0. 


( 2 . 12 ) 


From now on, the superscript indicates the species, and the subscript 
indicates slow or fast time-scale. The assemble average of Eq. (2.12) over 
the slow time-scale becomes 


a 

— (n s ) + 7 

at 


-t e e^ e 

n s v 8 + <n f v f > 


= 0 


(2.13) 


Eq. (2.13) subtracted from Eq. (2.12) becomes 


e-» e 


e-» e 


e-» e 


— (n f ) + 7 • (n s v f + n f v f + n f v s — <nfV f >) — 0, (2.14) 

at 

From Eq. (2.3), we obtain the lowest order component equation for 
the fast time-scale, as shown by 


3vi 


at 


(2.15) 


m. 


If we use the relationship in Eq. (2.15) to analyze the terms in, 
Eq. (2.14), the equation becomes 


e-» e 

7 • (n f v f ) 


— ( n f ) 

at 


e e 

k n f v f 

e 

uTi f 


k | e | Ef 

fk _i/2] 


W pe^ 

w m e w 

Uf 

lk d J 




(2.16a) 


and 
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v • (n f v s ) 

f k 1 


r wp e > 


rv s ] 

3 

ik d y 


l w j 


'^Te' 


— (n f ) 

at 

where k and w are the wave numbers and the frequency respectively, and k d and 
Upe are the characteristic wave numbers chosen to be the width of the 
current sheet and the electron plasma frequency respectively. The fast 
time-scale component energy density is defined by 


2 

E 

w f = — (2.16C) 

47m°T 
e e 

Corresponding to the slow time— scale fluid motion (i.e., collective 
motion) , 


. -* e . 

I V s I * V Ie 

with v Te representing electron thermal speed. 
Hence, in the case of 

V T e 

« 1/ (V^ = w/k) 

and 


(2.17) 


(2 . 18a) 


v Te 

« 

Eq. (2.14) 


f 

w _ - 1 / 2 
W f 

\ Wp e / 

could be simplified as, 


3n f _ e 

+ 7 • (n 8 v f ) 

at 


o. 


By combining Eqs.(2.15) and (2.19), 
approximate values: 


(2.18b) 


(2.19) 

we obtain the following 
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( 2 . 20 ) 


( 2 . 21 ) 


V w / 


Comparing Eqs. (2.17), (2.18b) (2.20) and (2.21) with the terms in Eq. 

(2.13) yields, 


e-» e 

n f v f 


n s v s 
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pe 


-1 / 2 

W, 


V u I 


“pe _l/2 
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V w 



V T e 1 

/ 

l|v e |J 


« 1, 


( 2 . 22 ) 


Thus, Eq. (2.13) becomes 
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-* e 


at 


(n 8 ) + v • (n 8 v s ) = 0. 


(2.23) 


For the ion, we may repeat our calculation by using Eqs. (2.12) through 
(2.22) together with the condition of Eq. (2.18b), which leads to a weak 
condition 




' p e 


\ w / 


Itl, 


« 


Vm i y 


w< 


.-1/2 


(2.24) 


' u pe^ 


Therefore, we obtained a similar expression of Eq. (2.19), as represented by 


at 


nj + 7 • (n s Vf ) = 0, 


From Eqs. (2.20) and (2.21), it could be shown, that 


1 

v f ~ 


and 


i 

''k ' 


f \ 

“pe 

-1/2 

'm e ' 

n f/n 8 - 
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IkJ 


V w J 


.m i; 


« 1, 


(2.25) 
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/ 
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U i; 




(2.26b) 


8 



Thus, we also found 


S -t i 

+ 7 • (n s v s ) - 0. 

at 


(2.27) 


By examining the electron and ion momentum equations according to Eq. 
(2.8) together with the principles of the order of magnitude approximation, 
we found the slow time-scale momentum equation for the electron as 


a -* e -* e -» e ® 

— V s + (v s • V) v s = — 

at 


-» -» e -» 

E s + - v 8 X B s 


m e n s 




-» e -♦ i -> -» e 

- v e i ( v s “ v s) + 9 + F p • 


(2.28) 


where F p is the pondermotive force due to the high frequency plasma andMHD 
interaction 


1 

2 


-t e -t « ty 

Fp 7 <(v f ) 2 > + 


-* e 


<V f X 7 X (0 e X B s )> 


m e c 


and 


V f = ^ e 

at 


(2.29) 


with being an arbitrary function. 

Similarly, we obtain the slow time— scale momentum equation 

for ion, 


9 -» i i x -» i “ e 

— V s + (V s • 7) V s = — 

at 


1 -i 


E s + - V 8 X B s 


mjn 8 




-> e i -» i 

+ ^ e i ( v s “ v s) + 9 + F p ' 


(2.30) 


where F p is the pondermotive force corresponding to the ion, 


1 e 

Fp = - — 7 <(Vf) 2 > - <Vf + 7 X (0 i x B s )>, 

* «■* 


mjC 


with 


9 


Vf *i- 

at 


(2.31) 


From (2.6), the fast time-scale field equation can be derived as, 


7 x B f = — 


1 3E f 4 ire 


c at 


-» e e-» e e-* e e-» e 

n s v f + n f v f + n f v s - <n f v f > 


-» i i-t i i -» i i-» i 

- n s v f - n f v c - n f v t + <n f v P > 


l f v s 


It also can be shown that 


L f V f 


LfVf 


(2.32) 


e-» e 


n f V 8 


U P € ^ 
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r i v h 

1 v s 1 

-» e 
n s v f 
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vk d y 


v T J 


« 1 


n s v f 


« 1 , 

n c 


so that Eq. (2.31) becomes 


v x B 


f - 


1 3E f 47ve _ e 
+ n s v f 


c at c 

Further, using the followwng expressions, 

e - * 

v f E f + - 

at 

and 


® -» e e -» 

- v f E f + 

m fl 


v f x B s 


m e c 


(2.33) 


(2.34) 


B* = - 


m e c 


d* t 


v x — - + v x (* e x B s ) , 


e at 

we derive the following equation: 

- 1 - 1 4 *e 2 _ 1 

7 x 7 x * e + — + — n s * e - — u Be x 

c i c i m e c z 


r& \ 


VB, 


(2.35) 


m e c 


7 X 7 X 


0 e x B s 


= 0 


(2 .36) 
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where 




^ ^ e -4 e 
V f 

8t 


“Be 


eB s 

m e c 


\ 




By defining 


P = n s (itti + m € ) 

(m iV s + m e V s ) 

u = 

(mi + m e ) 


Using Eqs (2.32) and (2.27), we obtained 
8 

p + 7 • (pU) = 0 

8t 


The LHS (left hand side) of Eqs (2.29) and (2 
as 


— (m a n 8 V s a ) + 7 • (m a n s V s a V s a ) , a = e, 
8t 

-» a -» a, -» cl 

where v s v s is the dyad of vector v s which is 

tensor. 

Since, 


n, 


pU u = 


e -» i -» e -» ) 

(m e v s + mjV s ) (m e v 8 + mjV, 


mi + m e 

= E (n s m a v 8 a v 8 a ) + m e n 8 [ (v 8 - v 8 ) 


in which, 


-t o 
Ds 


-4 e 

en 8 (v 8 



(2.37) 


(2.38) 

(2.39) 

31) can be written 
► • 

a second order 

) 

-4 e 4 i -* e 
+ V 8 (V s - V s ) ] , 

(2.40) 
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then, 


-t a a m e m s 

l n 8 m a v s v 8 = pu u + — — ]° ]° , 

e z n z 

s 


(2.41) 


By introducing the pressure tensor at the center of the mass 
system, we obtained 




P - n s [m e <(v - u) (v - u) 


-4 -4 


> + mj < (v - u) (v - u) 


f e 


fi (2.42) 


According to Eqs.(2.38) and (2.40), we have 
1 


-* e 

v s ~ 


en, 


1 0 -* 

Ds = u / 


-jo 


v s - u 3 


mi en s 


Incorporating these relationships into Eg. (2.42) yields 


-» e 


(2.43) 


P = m e n s <(v - v 8 ) (v - v s ) 


> + mjng < (v - Vi) (v - vj) 


f e 


f i 


1 + — 

ms ) 


m e o-4 o 4 e 4 i m e -* 0-10 
7 3 s3 s “ Ps + Ps T 3 s3 s 


(2.44) 


n„e‘ 


n s e‘ 


Multiplying Eqs. (2.28) and (2.30) by m e n s respectively and then adding 
them together results in 


/ 3u 


+ (u • 7) U 

at 


3 s x - 7P + P g + F . 


(2.45) 


-* i 


In deriving the above equation, we employed the relation F p << 


F p , so that 


p P = -- 


1 (m e ) 


m i y 


-» e 


/>v <(v f ) z > + p<Vf x 7 x (0 e x B s )>, (2.46) 


n^c 


In addition, we employed the general Ohm's law similar to those given by 
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Krall and Trivelpiece (1973) by combining Eqs. (2.28) and (2.30) and 
multiplying the terms of e/m e and -e/m e respectively In case we neglected 
the second order terms together with the pressure gradient, gravitational 
force and pondermotive force, we get 


-* 0 -* -» -» 

*>3s ■ E 8 + U X B s/ C 

where 

t) = m e K ei /(n s e 2 ) 

From Eg. (2.6), the slow time-scale current becomes 


(2.47) 


(2.48) 


“t 

Us 


_> e i e-» e i-» i 

en s ( v 3 “ v 8 ) + e <n f v f - n f v f > 


(2.49) 

According to Eq. (2.26) , the second term on the RHS (right hand side) of Eq. 
(2.49) can be delected, thus 


-1 _ 7° 

3 s 3 s 

-t e 

= en s (v 8 - 

-t i 

V s ) 

(2.50) 

Therefore, 

-4 




1 3E s 4 7t 

. 0 47V 0 


7 X B s 

S 

i 0 i u 

3 S = — 3 S 

(2.51) 

c at c 

c 



We dropped the displacement current because << «p e ■ 

Finally, we derived a set of global coupling MHD equations 
with pondermotive force, 


p + 7 • ( pU) = 0 

at 


. au 

p| — + (u • V) u 

at 


^ -» ■» -* ■i, 

= - j x B - 7p + P g + Fj 
c 


(dh 


7 X 7 X 


lat } 


1 a 3 * 1 4 ire 2 dt> 


c at 3 c 2 m e m 

( w Be ® 2 ^ 

'B 1 

r x 

— 

l c 2 at 2 

IbJ 
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m e c 


7 X 7 X ( 0 X B) 


(2.52) 


(2.53) 


= 0 

(2.54) 



with 


e ^ 1 m e 

F p = p <4> x v x (0xB)>- — — p v < 


n^e 


2 m; 


r a^ 

^at, 


(2.55) 


Ohm's law, 

1)3 = E + U X B/C 
Maxwell equation, 

* 4 * i 

7 x B = — J 

c 


7 X E 


1 3B 

c at 


(2.56) 


(2.57) 


7 • B = 0 

/ 

In the case of w pe » u> Be ; « » w Be and * I I® (i-e- / the hi 9 h frequency 
amplitude parallel to the magnetic field), the Eq. (2.54) and (2.55) are 
reduced to 


_ 1 a 1 2 v f 1 4*e 2 _ 

7 X 7 X V f + + — — pV f = 0 

c 2 at 2 c 2 m e mj 

1 m e _ 

F p p 7 <(v f ) 2 > 

2 m i 


(2.58) 


(2.59) 


with energy equation, this set of equations form a closed set of governing 
equations for the described problem. 

2.3. Resistive Instability by Soliton 
It has been recognized that one way to convert magnetic energy to kinetic 
energy is through the resistive instabilities. Although the foundation of 
this study has been made by Furth, Killeen and Rosenbluth (1963 thereafter 
it will be referred to as FKR) , Priest (1986) , the analysis including the 
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pondermotive force is still lacking. Thus, we shall present a detailed 
analysis for the resistive instability by soliton for a transverse plasma 
wave in a slab geometry. 

For convenience, we shall perform an incompressible analysis, which 
can be justified in the most physical region (Furth, Killeen and 
Rosenbluth , 1963; Shivamoggi, 1985) . The motion of plasma is along the x- 
direction and the unperturbed field is B oy in perpendicular to the x- 
direction, and assuming that the gravitational effect is negligible, we 
have the governing equation as follows; 


3 P 

+ 7 • ( pU) = 0 

at 


3u 

— + (u • v)u 

l at 


(3.1) 


1 m e 2 

= — (vxB) xB-vp- pV<V f > 


1 

4tt 


2 mj 


3B 


vc‘ 


7 2 B 


1 4ite‘ 


= 7 X (U X B) + 

at 4* 

1 3 2 V f 

7 X 7 X V f + + V f - 0, 

c 2 at 2 c 2 m e m t 
1 3B 

7 x E = , 

c at 


■* 4 * , 

7 x B = — 3 , 

c 

7 • B = 0 . 


(3.2) 

(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 


Let us examine the steady state of this set of governing equations (i.e, 
Eqs . (3.1) through (3.2)) by setting 


u 0 = 0 


p 0 = p + p*(X) 


(3.8) 
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B 0 = B oy (x) y = B tan h 


'x ' 


(3.9) 




Within the current sheet ( i . e . , | x | - | x/L, | « 1) , the magnetic field can be 
approximated by 

B 0 = B — — y < 3 - 10 > 

Ls 

where L s represents the characteristic length of the current sheet. 
Incorporating Eg. (3.8) into Eg. (3.2) , we obtain 


1 c s 1 m e -> 2 

— — (7 X B Q ) X B 0 - — 7p' 7 <V f0 > 

4 it p P 2 m i 

where c s is the lowest order sound speed. 


= 0, 


'^p^ 


c s = 


dp 


Using Egs. (3.10) and (3.11), we obtain 
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p 

m e 

r>2 

2 B 

X 1 


= " ? 

[2 


<Vfo> + — 
8tt 

VL S ; 



(3.11) 


(3.12) 


(3.13) 


B 2 


\ 2 


In the case of, 
1 

8 tv 

we have, 


lL, 


1 _ 2 

« - P <Vfo> 

2 mi 


1 9 m e -,2 


<V f „>• 


(3.14) 


(3.15) 


c' 2 mi 

S 

As usual, the lowest order of fast oscillation speed of the electron can be 
expressed by 


v f o = “ 


1 

2 


v 0 (r,t) e 1M t + 


(3.16) 
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With w 0 being the lowest order fast oscillation frequency. From Eq. 
(3. 16), it can be shown that 
„ 1 
2 


<v f0 > = - |V 0 


(3.17) 


i a 

lnv, 


« 0 5t 


« 1 


(3.18) 


and by neglecting the second order time derivatives, Eq. (3.4) 
becomes 


2iu, 


3v, 


at 


2 2 " > 

- C 2 V 2 V - 


' 2 4Tve 2 Po ) 


m B nt I 


v 0 = 0. 


(3.19) 


For the transverse plasma wave, 


2 2 


4fle 2 p 


pe 


(3.20) 


m e mi 


Employing Eqs. (3.15) and (3.17), Eq. (3.19) becomes 


3v, 


at 


v z v 0 + 


* p e 


w pe m e -* . , -> 

— |v 0 | 2 V 0 = 0 , 

8c z m ; 


(3.21) 


For convenience, let us take the complex conjugate of the Eq. (3.19) which 
leads to a standard nonlinear Schrodinger equation. 


3v 0 

i = 

at 

and where 
v„ = v, 


i a z v 


0 1 i 9 

6 v n v, 


(3.22) 


2 3X‘ 


(x : t) y, X = - -jr X 


S = 


“pe 

8c 2 


U , 


u = 


c 

m e 

mi 


pe 


(3.23) 


From Eq. (3.23) it is clearly indicated that these nonlinear wave 
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interactions can be considered as wave packets and its wave function is 
proportional to v 0 . These wave packets can produce self-generated 
potential, B| v 0 1 2 . Since 6 is positive, this self -generated potential has 
attractive and repulsive characteristics. In other words, these 
nonlinear wave packets can be self-centered and form a stable structure 
called "soliton". For the steady state, the solution for this soliton can 
be expressed by 


v n = v n sech 


■J5 v 0 x 




(3.24a) 


6 

i> = — t + 0 O 
2 


(3.24b) 


By incorporating Eqs. (3.24a) and (3.17) into Eq. (3.15), we 
obtain 

i . ( x ) 

, (3.25) 


Po ( x ) = " 


P 0 

— u(v 0 ) Sech 


C ‘ 4 

s 




the width of the soliton being 
8 ( c 

S u- 


f o = 


W C s ' 


V w 


(3.26) 


pe 


o 


Hence, the width of the current sheet can be estimated from 
Eq. (3.14) , such as, 


AX = € « 


rrr 

4 

o 

. 


* V, J 


(3.27) 


with v A being the Alfven speed. 

Now, let us turn our attention to the perturbed state. By 

using Eqs. (3.1) through (3.3), we obtained. 
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3u 1 .» ^ -* -* 

p 0 = _ 7 p + (7 X B) x B 0 + (7 X BJ X B 


at 


47V 


47V 


ap7|v 0 


(3.28) 


3B _ _ 1 , 

= 7 X (U X B 0 ) T)C Z 7 X (7 X B) 

at 47V 


dp 


= - (U • 7 ) p 0 

at 


(3.29) 

(3.30) 

(3.31) 


7*11 = 0, 7 • B = 0 

By taking all the perturbed quantities, such as 

A(x,y,t) = A(x)e 7t e l ky (3.32) 

and substituting these perturbed quantities Eq. (3.32) and Eqs. 
(3.28) through (3.30) together with Eq. (3.31) to eliminate B y , 
and u y , we obtained a set of differential equations, 


d ( , du x ) 
7 P o 


dx 


dx 


- k' 


1 dP„ d| v. 


7 Po u 


4 dx dx 


f)C' 


B oy (x) 


u. 


ik 

4tv 


~ B 0 y t x ) 7 


47V 


iC 


1 d 2 B 


1 - 


oy 


B oy dx' 


By r 


1 B x = ikB oy (x)u x + 


nc 2 ( d 2 


47V 


-k 2 


ldx‘ 


By# 


(3.33) 

(3.34) 


Introducing the following dimensionless quantities, 

B x 

B oy (x) = B F (X) , i = — , Po = P e o( x )» 


{ 1 d 

G 0 = - - U — |v, 

l 4 dx 


B 

' 1 


- P n t h 
0 
p 


(3.35a) 


(3.35b) 


X = L S X , a = kL s , S - 


(3.35c) 
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To = 


47TL, 


t»C‘ 


Tu = 4 7tL 


^7 

s — ' 
B 


( 3 . 35d) 


and incorporating them into Eqs. (3.33) and (3.34) yields. 


d 2 * 

_ 2 , 

- 3 — = a + 

dx 2 


d ( du x ' 


1 + 




a 


- i k t r F u x , 


(3.36) 


dx 


'0 - 

dx 


= a 2 u. 


So + 


1 2 T 2 


G 0 + 


7 t , 


r i 1 

( 1 d 2 F) 

iff cl S r “ 

l 7 T r dx 2 J 

, kr R J 


(3.37) 


In order to seek solutions for this set of governing equations, we 
shall divide them into two regions; the region outside the current sheet and 
the region inside the current sheet. 

i) Region outside current sheet. 

In the case outside the current sheet, the -n 0 in this region; thus, 
s 2 00, under this circumstance, Eq. (3.37) becomes 


f 1 


y 2 T 2 


G 0 + 


y t 


u 


-1 

x - * 

kr b 


( 1 d 2 F 1 

F 


*yt r dx 2 


rv 

Combining this experssion with Eq. (3.36), we get. 


ikT 0 u 


R U X 


dx 2 


1 d 2 F'\ 

a 2 + - — 

F dx 2 J 


F 7t r 

On the other hand, Eq. (3.36) can be rewritten as 


y r 


1 ( d2 * , ' 
_ - - a 2 <J> 

,dx 2 


« 1 - ik - u x 

y 


(3.38) 


(3.39) 


It is also true that 7 t r -» <» outside the current sheet, thus 
Eq. (3.39) reduces to 
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(3.40) 


l 2 ik - Fu x 
7 


(3.41) 


Combining Eqs. (3.39) and (3.40), we obtained 

d 2 * ( 1 d 2 F G ^ 

— — — - ^ a 2 + — _ + — = 0 , 

dx 2 l F dx 2 F 2 J 

Concerning the transverse plasma wave, G 0 is negligible outside the current 
sheet because the incompressible characteristics are held outside the 
current sheet; therefore Eq. (3.41) becomes 


d 2 * 


dx' 


- * 


1 d z F 


a £ + — 


F dx' 


= 0, 


The solution for this equation is 


(3.42) 


*+ “ 


. -ax 


.ax 


tanh x 


1 + 


1 - 


a 


tanh x 


a 


x > 0 


x < 0 


The jump condition across the current sheet is 


A' = 
e 


1 

d/ 

d^ 

’ 

= 2 

i 1 

— - a 

rf(0) 

^dx 

+ o 

-0 J 


iCL J 


ii) Region inside the current sheet 


(3.43) 


d 2 u„ de 0 

Since, — — >> — =: inside the current sheet, then Eq. (3.37) 
dx 2 dx 

becomes 


0 , 


ct 2 s 2 


7 T . 


d 2 u 


- 


a U 


dx' 


1 G, 


— + F 


7T R F 


( d 2 ^ 


_ 


a </> 


dx' 


d 2 F 


- + - + — — , 

F dx z 

and Eq. (3.36) can be rewritten as 


(3.44) 
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1 (d z * 


7T R ldx‘ ) ^T R 

where 

u = - i k T r u x 

Using Eqs. (3.8), (3.10), (3.25) and (3.35), we find, 

d 2 F 

F - x , — — = 0 


(3.45) 


1 d r f 0 ) 2 X p u. ( 0 ) d 

G 0 = - - u — v 0 sech 2 — - — v o - sech 

4 dx 1 J e oJpL c „ 4 ^ ^ 


u 2 fv7 2 o) 2 [ d f X 

16 Ic J l J H L dx l U 0 


J s 1 — 


= 7V il { 


By noticing the form of Eqs. (3.8) and (3.25), we may choose 
9 0 Z 1; thus Eqs. (3.44) and (3.45) become 

1 ( d 2 u ) _ (d 2 * „ ) 


1 

“V T 

d 2 u 

a 2 S2 

7 1 R 

\dX2 

1 

'd 2 / 

idx 2 

\ 

- a 2 ^ 

/ 

with 



£ = 

do 



a 2 U = (1 + £) X — — " a- 2 S ~ ?^( aT R) X, 

j-- 2 


U X . 


d n = TV u i — 
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In the case of highly conducting plasma medium, then we can apply r R » 
1 and s » 1; in addition if the characteristic scale length "x" is much 
smaller inside the current sheet in comparison with any other 
characteristic length, then we may apply the theory of resistive ordering 
suggested by Dobrott et . al . , (1977) , thus; 


7 T R ~ 5 ■ 1 , X- 5, ^ - 1; 

From the RHS of Eq. (3.48), we note, 
* - 8 2 u, 
which leads to 
u ~ 5 “ 2 . 

From Eq. (3.49), we find 

1 d 2 u _ 

(7t r ) — - U(7t r ) x 
dx 2 


(3.51a) 


(3.51b) 


a 2 S 2 

yielding 

s 2 ~ 5 ~ 5 . (3.51C) 

From these relationships, we observed that both functions of i and u are 
proportional to 5 ; thus we may expand these two functions as a series of 
functions respectively, 


^=^0 + ^l + ^2 +, ‘* 
U = Ug + U J + Uj t • « • 

In Eq. (3.48), 


(3.52) 


1 

( II II 

. ' 1 


f \ 

— 

^ — CL yp — a i + • • • 

3S 

+ ^ l + * * * 

^R 


j 


i / 


+ 



Ug + U 1 + 



7Tg ( ) 

where superscript primes represent derivatives with respect to 
Using the order of approximation, we note that 


x. 
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t‘ - , *; - > 

0 1 

. - 2 


5 


2 _ 


5 


u 0 - a # u x - a 

which leads to 

,■ , o < 3 - 53 > 

0 

This condition gives constant approximation, 

<l> Q = constant (3.54) 

Based on the expansion given in Eq. (3.52), we obtained the following 
successive order of approxmation governing equations: the first order 
equation for / is, 


1 1 

- * 0 + 

-»7 R 1 IT, 


u 0 x 


(3.55) 


Similarly, we obtained as the first order equation for u from 
(3.48) 

7 T 


a 2 S 2 0 


u" = - s (*vt r ) x + (i + n x i . 


(3.56) 


Combining Eqs. (3.55) and (3.56), we have 

y t 


a 2 s 2 0 


u" - (1 + £) X 2 u 0 = (7T R ) x^ 


0 • 


(3.57) 


From Eq. (3.55), the jump condition within the current sheet can 
be computed as, 


ft 00 


A' = 

1 


1 P 


dx = — 

1 


( 7 T R + X U 0 ) dX 


Assume that 

u 0 (x) = *h(x) : x = Ex 

where 


(3.58) 


(3.59) 
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(3.60) 


7T, 


l A = and k - 7t r /L 


a 2 S 2 


using Eqs. (3.58) - (3.60), Eq. (3.57) becomes 
h" (X) - (1 + n x 2 h(x) = * 0 x ; 
Correspondingly, Eq. (3.58) becomes 


Ai - — (7T R ) 5/4 (as) 

K 


- 1/2 


,co 


Uo + xh(x)) dx 


Introducing the Fourier transformation, 


h = 


00 


,-lZX 


h (z) dz. 


and its inverse transformation, 


1 P® 


h(z) 


2 TT J _ 


. i z x 


h(x) dx. 


00 


Then Eq. (3.61) becomes, 


h " ( z ) - (1 + ?) h(z) z 2 = - / 0 5 1 ( z ) , 

i 


(3.61) 


(3.62) 


(3.63a) 


(3.63b) 


(3.64) 


where 5(z) is the Dirac Delta function; hence Eq. (3.64) reduces 
to 


h“ (z) - (1 + 5) h(z)z 2 = 0, for z f 0. 


(3 . 65) 


The solution of the above equation could be represented by the MacDonald 
function (Krall and Trielpiece, 1973) . 


* + 


h (z) = A^Jz K 1/4 


(1 + £) 1/2 z 2 /2 


, z > 0, 


(3.66a) 


h'(z) = -iA^Jz K 1/4 (1 + %) UZ z 2 /2 


, Z < 0. 


(3.66b) 


The appearance of the solution (Eq. (3.66b)) is due to the con- 
sideration of the property of an odd symmetry for the RHS of Eq. 
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(3.64) . 

For the region z = 0, Eq. (3.64) becomes 


» II ^ o 

h = — 8' (z) 
i 

By integration of both sides, we have 


(<j, \ 


\i ) 


2A = 


H + (Z) 


(3.67a) 


Z->0 


with 


H + (z) = -nIz K 1/4 


(i + n 1/2 — 


(3.67b) 


Since 


p® 




- 00 


r>°° 


8 (z) dz 


- 00 
P®. 


dxe" 1 zx = 2 tt^ 0 8 (z) 


xh(x)dx = 


h(z) dz 


- 00 


( a p® 


1 9z J _ 


le 


1ZX dx 


z = 0 
p® . 


= 2iri 


h(z) 8 ' (z) dz 


= - 2 7Tvf o 5 ( Z ) 


- 2 7V1 


P® . 

h(z) 8 (z) dz , 


z=0 


+ oo 


P® 


J 0 


h» (z) 8 (z) dz = 2 

00 

We finally derived 


d 

— H + (z) 
(7r R ) 5/4 dz 

A ' = - 2tt 

1 (as) 


h' (z) 8 ( z ) dz . 


z=0 + 


1/2 H + (Z) 

using the characteristics of the MacDonald function, 
d ' 


dz 


z 1/4 K 1/4 (z) 


= - Z 3/4 K 3 / 4 (Z) 


(3.68) 


and the asymptotic value for |z| << 1, 
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TV 2 1 

Z K K k (Z) : ' 

2 sin v tv r ( 1— k ) 

with T being the gama function, Eq. (3.68) becomes 
(7t r ) 5/4 r(3/4) 


(3.69) 


a- = 4 7V (i + n 

i (as) 1/2 T(l/4) 


1/4 


(3.70) 


By equating the jump conditions inside and outside the current sheet, we 
find the growth rate dispersion relation from Eqs. (3.43) and (3.70), 


(7t r ) 5/4 r(3/4) 1 1 

47V — a + u 1/4 

(as) 1/2 r(l/4) 2 U 


(3.71) 


This is the key equation for determining the resistive instablility . 

For the illustration, we take a coronal structure of current sheet, 
the typical plasma and field parameters are: 


10 km. 


B = 100 Gauss, 


n e = 10 9 /cm 3 , T e = 10 6 ok. 


(3.72) 


Based on these fundamental values, we obtained Alfven speed (v A ) 
to be 8.4 x 10 7 cm/s and the sound speed (c s ) to be 2.3 
x 10 7 cm/s. If we take v° ~ 3.10 7 cm/s, from Eq. (3.50b), we 

find that, 


d 2 = 2.2 x 10 4 . 
0 


(3.73) 


Correspondingly, the resistive parameters are; t r ~ 1.3 x 10 6 s, 
Th = Ls /v a - 1.2 x lO' 2 s, s = t r /t„ = 10 8 . Therefore, Eq. 

(3.71) gives 

-y . ^1/2 ( 3 - 74 ) 
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which says that the growth rate is proportional to the square root of 
resistivity, such as 


7 


1 a 1/2 s 1/2 r(l/4) 1 


d 1 ' 2 



\ 

a 


r(3/4) 


t r 


3 • 10 -3 (S - 1 ) . 


(3.75) 


This result is based on constant solution (i.e. regular homogeneous 
field configuration) . When the reconnection process begins, the constant 
approximation no longer becomes valid; thus, the other kind of 
instability may appear. We shall examine this situation in the next 

section. 


2.4 Eruptive Local Instability 

When the constant approximation becomes involved, we should re- 
examine the set of coupling equations (i.e. Eqs. (3.48) and (3.49)). By 
introducing the Fourier transformation as suggested by FKR, such as, 


t = 


;(fl)e 


“ i8x do 


(4.1) 


u = 


i0x do 


u( 0) e' 

CO 

into Eqs. (3.48) and (3.49), we find that 


(4.2) 


7 T i 


a 2 S 2 


(0 2 + a 2 ) U(0) « (1 + £) 


d 2 u(0) dvf(0) 


+ 17 T i 


d0 


d0 


(4.3) 


1 1 du(0) 

(0 2 + a 2 ) ^(0) = - /( e ) + i 

7 T p dO 


7 T i 


(4.4) 


d\ff 

Using Eq. (4.4) to eliminate — in Eq. (4.3), we obtain 

do 
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0 


p 


d 2 u( 0) d 

*f 

do 2 do 


r e 2 + a 2 ' 
i O 2 + a 2 + P; 


du( 0) 

do 


u(0) 


la - 


+ 1 


(4.5) 


with 

(4.6) 

P = 7t r 

It is worth noting that the Eq. (4.5) is identical to the Eq. (5) of FKR by 
setting 5 = 0, which implies that 7 t r * ® or d 0 ■> 0. Physically, this 
corresponds to a very thin current sheet and a high resistive region. 
Under the present case, we have 


— = a « 1 ? 
P 


(4.7) 


hence Eq.(4.5) can be written as 


d 2 u(0 2 ) d 
5 — + 


f 0 ! + <r du ( 0 i ) ) 


dO' 


do 


+ 1 do 1 


A(9 2 + cJufOj) 


(4.8) 


with 

0l = p-i/2 d and A = P 3 / a 2 s 2 ( 4 - 9 ) 

Eq. (4.8) is the equation to be used to determine the solution u(0) and 
subsequently ^ will also be determined by Eq. (4.4) . Finally, the growth 
rate for the non-constant-; -solution could be found. 

in order to seek the solution for Eq. (4.8) , we introduce the following 


parameters : 


0 


<r 1/2 


0 

a 


(4.10) 


- LJ s 

Since a = — << 1/ thus, 

L y - 

(14.8) becomes x -> ± ®, 


>> 1. In this limiting case, Eq. 
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A 



d 2 u 

d 1 

A 

du 


— A •— "H “ A “ 

(1 + « 2 ) 

— 

c 

dO 2 

d0 


do 


For £ > 1 , we have 


A A 

A<r( 0 2 + 1 ) U = 0 


(4.11) 


£ 9 2 * , 

— » : 8 2 

<r <r 

Eq. (4.11) gives 

A A 

% d 2 u . du 

+ 20 A o - ( 0 2 + 1) u = 0 (4.12) 

<r d0 2 dO 

The above differential equation can be simplified by the following 
transformation 


<r * 


u ( O ) = exp 


2$ 


w ( 0 ) , 


so that, Eq. (4.12) reduces to 


d 2 w 


A - (B 0 + 0") w = 0 


do 


where 


5 i' 

A 0 = — (£A + I)" 1 


f(l + A<r) 


B 0 = 


<r(£A + 1) 

Eq. (4.14) can be simplified further by introducing 

-J2 » 


y - 


1/2 


9 / 


(4.13) 


(4.14) 


(4.15) 


(4.16) 


so that the final form of Eq. (4.14) is 


l 
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(4.17) 


d 2 w 

+ (g + 1/2 - 1/4 y 2 )w = o 

dy 2 


where , 


B, 


g + - “ T77 

2 2A 1 ' 2 


(4.18) 


It is well known that Eq. (4.17) has the solution of a parabolic cylindrical 
function thus, 


w = D g (y) 

when g = -1, D g (y) 


(4.19) 


satisfies the following boundary condition; 


D„ (y) -» 0, as | y | » 1. 


Thus, 


7t 

D_! (y) = N2 


exp 


'1 

- y 

1 4 ) 


i - * (y/^2) 


(4.20) 


( 4 . 21 ) 


Where B (y) is the probability integral, its asympotic value becomes 


1 2 n - 1 (-D ra r(m+l/2) 

1 - HV/-&) m - e ' Y 2 


m = 0 n\in+l/2 


\ 2 ) 


j y | » 1 , (4.22) 


.2 m + 1 


Returning to Eq. (4.12) for g — -1 implies 
B„ 


, 1/2 


= l . 


(4.23) 


or (1 + A<r) 2 — (1 + £A) 


(4.24) 


Using Eq. (4.24) together with Eqs. (4.7), (4.9) and (3.50), we 


we found the growth rate for non-constant -4 to be 
2 


7 = 


- 2 


(4.25) 


la' 
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2 . 5 Numerical Results 


In order to illustrate our theoretical results, we have chosen a 
coronal structure of a current sheet plasma for testing. The typical 
physical parameters for this kind of plasma and field are given in Eqs. 
(3.72) and (3.74), which leads to 7 = 3 X 10“ 3 s' 1 for the constant case. 
As a result we use the very same values for Eq. (4.25), yielding 7 2 1.2 X 10 s 
s -i for the non-constant -* case. It is immediately reconized that the 
growth rate between these two kinds of instability has a 10 7 order of 
magnitude. Thus we can conclude that the constant case leads to a 
resistive instability of slow time-scale phenomena and the non-constant 
case leads to a much faster time-scale phenomena called "eruptive 
instability". Physically, we may claim that when a magnetized plasma 
begins its dissipation in an orderly fashion through resistive 
instability, the inbedded field is rather uniform. As soon as the field 
configuration is distorted, the energy release can be more violent because 

of eruptive instability. 

2.6 Concluding Remarks 

In this paper we have presented a theoretical result concerning the 
transition of resistive instability to eruptive instability with 
pondermotive force with high-frequency plasma ossicillations. It shows 
that the constant -ix (i.e. uniform field configuration) leads to resistive 
instability and that the non-constant -i case leads to eruptive 
instability. The close form solutions based on plasma and field 

parameters are given for these two cases. 

A numerical example is also given showing that the growth rate for 
these two cases can be different in 10 7 orders of magnitude for the solar 
corona structure of the current sheet. 
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3 ANALYSIS OP INSTRUMENTAL POLARIZATION OF IMAGING OPTICS 
3.1 Int roduct ion 

The objective of the SAMEX magnetograph's optical 
system is to accurately measure the polarization state of 
sunlight in a narrow spectral bandwidth over the field of view of 
an active region to make an accurate determination of the 
magnetic fields in that region. Our design goal is to measure 
magnetic fields to an accuracy of one part in 10 4 . To achieve 
this accuracy requires a polarimetric accuracy of 10 in 
determining the polarization components of the light as a 
fraction of the total intensity. This requirement means that the 
instrumental polarization of the optics must be reduced to levels 

below 10" 5 . 

All optical elements introduce some polarization 
change, especially when used off axis. Combinations of mirrors 
and antireflection-coated lenses can display a full range of 
polarization behavior: linear and circular polarization and 

linear and circular retardance. Of particular concern in the 
development of highly accurate polarimeters is polarization 
rotation which causes linearly polarized light to leak through 

subsequent of crossed polarizers. 

Because of this instrumental polarization, the SAMEX 

foreoptics (the optical elements in front of the polarizer - the 
Cassegrain telescope and the relay lenses) must be considered as 
a weak polarizer in front of the polarimeter section. This 
instrumental polarization in the foreoptics changes the 
polarization state of the sunlight incident on the polarimeter 
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and thus introduces errors in the measurement of the . solar 
maanetic field. 

If possible, it would be preferable to place the 
polarimeter before the imacjina oDtics, to locate polarizers and 
retarders in front of the Cassegrain telescope. Then the 
determination of the polarization state of liaht would 
unaffected by the instrumental polarization of the foreootics. 

This desian would have the additional advantaae of allowing the 
liaht to pass throuah the polarimeter over a smaller ranae 
angles of incidence. However, such a design is impractical for a 
system with a 30 cm aperture - high guality polarizers and 
retarders have much smaller apertures, on the order of a few 
centimeters. Thus, it is necessary to use collecting optics to 
collect the 30 cm aperture of light and focus it through small 
polarizing elements in the polarimeter section. As the collected 
light is passed throuah smaller apertures, the angular spectrum 
of the light increases in a relationship governed by the Lagrange 
invariant. By reducing the beam from 30 cm to 2 cm through the 
tunable filter and polarimeter, the range of angles-of incidence 
is increased by a factor of 15, from 5 min of arc to 75 min of 
arc, thereby increasing the anale of incidence effects in the 
lenses just in front of the polarimeter and in the polarimeter 
itself. Conseauently we anticipate there will be some 
instrumental polarization due to the foreoptics in the light 
entering the polarimeter. 

The goal of the analyses outlined in this section is to 
precisely characterize the extent of this instrumental 
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polarization, and to design the optics and coatings in unison to 
minimize this spurious polarization introduced by the 
foreoptics. We will show calculations of the instrumental 
polarization of ordinary foreoptics compared with the 
polarization performance of special ultra-low polarization 
optical coatinas desioned for this application. 

The instrumental polarization analysis uses a program 
tha incorporates the theory of polarization into the standard 
geometrical optics and lens design codes (Chipman, 1987). By 
including the polarization of the optical elements in the *=xrst- 
and third-order desicn process, the effects of coatings on curved 
substrates can be treated. For each ray, a polarization matrix _ 
(in the Jones matrix formulation) is calculated for the ray at 
each optical interface. These matrices are multiplied together 
to calculate the polarization matrix for that ray from object 
space to image space. An analytic function for the rays provides 
the polarization behavior as a function of the exit pupil, object 
height, and wavelength. This technigue represents a quantum jump 
improvement in the practical design of foreoptics in front of a 
polarimeter and will enable the degree of residual instrumental 
polarization to be reduced to < 10" 5 in the SAMEX magnetograph 

system. 


3.2 Instrumental Polarization of Standard Cassegrain Telescopes 
To establish the need for a detailed polarization 
analysis, we estimate the polarization effects associated with a 
Cassegrain telescope with aluminum thin film coatings. The 
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electric vector of light incident on a surface of the Cassegrain 
can be decomposed into two parts, the components of the vector 
vibrating parallel (p) and perpendicular (s) to the plane of 
incidence. These components have different reflections as a 
function of the angles of incidence. Figure 23. a shows the 
reflectance for s and p polarized light from an aluminum thin 
film coating with complex index of refraction n = 0.7 - 7.0 i. 
Figure 23. b shows the phase change on reflection for the s and p 
components. The difference in s and p reflectance causes a weak 
linear polarization aligned perpendicular to the plane of 
incidence to be associated with reflection from all metallic 
interfaces. Moreover, the differences in phase change cause a _ 
weak linear retardance to be associated with reflection from 

mirrors . 

The percent reflections from the aluminum surface for 
the surface parallel (s) and surface perpendicular (p) rays are 
approximately given by 

R = 0.94 + 0 . 02i 2 , 


and 

Rp = 0.94 - 0.02i 2 , 

where the reflectance rate of change per unit angle of incidence 
squared (i 2 ) is approximated from Figure 23. a. Therefore, the 
induced linear polarization is estimated to be 
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Fiaure 23. Linear polarization and retardance effects for a 
Cassegrain telescope. The percent reflectance (a) and phase 
chanoe (b) on reflection from an aluminum thin-filra coating is 

for bSSh the s and p components of the incident light wave. 
The wavelength of the light is 5250 A and the complex refractive 

index of the aluminum coating is n - 0.7 - '• Tn ® . . . 

differences in the s and p reflectance cause linear polarization 
aligned with the incident plane. The differences m the s and p 
phase cause a linear retardance. The differences are small for 
small angles of incidence, but they are. not negligible for the 
SAMEX magnetograph design. 
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At the edge o£ the primary mirror of aperture D and focal length 
£1, a ray is reflected through a total angle of tan(0/fl) ~ D/,£1 
„ l/£, where £ is the f-ratio. Then for a Cassegrain telescope 
the maximum angle of incidence i„ is given approximately by one 
half the reciprocal of the f-ratio. for an f-ratio of 5, 

angle of incidence is i m *l/10. Hence the induced linear 

a nf Q ♦ 02 a 2 x 10~ 4 . For a linearly 
polarization is on the order of 2 x iU 

polarized ray this represents a rotation- of its plane of 
polarization over a segment of the mirror, and this rotation ^ 
Introduces errors in the deduced polarization state. In Figure 
24 the angle of incidence versus the pupil coordinate for an 
■Illustrative' Cassegrain telescope is shown for an on-axis and 
off-axis ray. The point to note is that, for the off axis ray, 
the average angle of incidence Is not zero, and hence there is a 
net polarization effect associated with the off-axis rays. 
Furthermore, even on-axls, we average over the square of the 
incident angle which gives a net polarization contribution, as 

shall be discussed below. 

For an on-axis beam incident on either of the aluminum 
coated Cassegrain telescope mirrors, the magnitude 
orientation of the linear polarization and linear retardance 
associated, with the mirrors have the forms shown in Figure 25.a 
(polarization) and Figure 25. b (retardance) as a function of the 
pupil coordinates. In Figure 25. a, the linear polarization is 
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ON AXIS, MIRROR 1 

ON AXIS, MIRROR 2 

OFF AXIS, MIRROR 1 

OFF AXIS, MIRROR 2 
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Figure 24 show 0 the%ariat ions™^ tL angle 

Cassegrain, mirrors. These cur coordinate for on-axis 

° f d i o1f d !!lis rays^ incident on both the primary and secondary 
mirrors of a Cassegrain system. ^“^fp^ry tn5 secondary 

•zero average angle of incidence at both the prim £ with 

mirrors. This implies a net li. " e « shoun heC e Is tor a 

the telescope *“ larger Afield if view than the one chosen in the 
SAMEX^ design^ th It ilTcSoUn to provide a vivid example of the 
off-axis problem. 
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oat terns * are* quadratic fSnltlons of the pupil coordinate and are 
termed -polarization defocus. " 
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zero in the center of the mirror where the beam is at normal 
incidence. The magnitude of the linear polarization 
(polarizance ) increases quadratically with distance from the 
center of the mirror. The linear polarization is oriented 
tangentially. This polarization pattern (polarization 
aberration) is called linear polarization defocus (or quadpol). 
Figures 26. a and 26. b show the effect of large amounts of linear 
polarization defocus (about 40 times more than the conventional 
Cassegrain) on beams of uniform lef t-circularly polarized light 
(a) and uniform, vert ical-1 inearly polarized light (b). The 
linear polarization associated with the telescope mirrors changes 
the polarization state of the light causing spatial variations of 
intensity and polarization across the beam. If the polarization 
state of the light is now measured with a polarimeter, a 
polarization state different from that incident on the mirrors is 
obtained. Despite the symmetry associated with the resulting 
transmitted polarization patterns (such as in Figure 26), the 
polarization variations do not cancel (due to averaging over the 
incident angle- squared) . This is best understood by considering 
the transmitted light as being a superposition of two 
polarization states: the incident state yields the correct 
polarization measurement; the light in the orthogonal state 
constitutes the error signal introduced by the instrumental 
polarization . 

For example, in the case of Figure 26. b, the light in 
the orthogonal state (horizontal linear polarization) has the 
form across the pupil shown in Figure 27. Although the phases of 
this horizontal component are 180° out of phase in the four 
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Figure 26. Linear golarization , 

telescope. This figure lllus ^5 t *t r rirfor a beam of (a) uniform 
effects for a Ca * se 9***? £ r a ® d of (b) uniform vertical Unea * 
left circular s a „ d arrows at the top, bottom, let t, 

polarization. The ellipses same relative position of 

ss sffiiM: 

Cassegrain** ‘STSSS-?* 

‘ phase of the light where one cycl i^ o£ polarizat ion which 
effect introduces the orthogo roduC ed by the instrumental 

llllZl SoTiTS: measured. nts^ made 3th the ^agnetograph. 
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* * signal of the magnetograph. 
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quadrants (resembling astigmatism), all this light will pass 
through a horizontal linear polarizer and contribute to the 
polarization error signal. The phase differences do not cause 
cancellation of the polarization aberration; instead they affect 
the polarization' accuracy and change the structure of the 

diffraction pattern (Kubota and Inoue, 1959). 

For a metallic mirror, the effects of the linear 

retardance are orders of magnitude larger than the effects 
weak linear polarization. Figure 25.b shows the form of the 
linear retardance associated with an on-axis beam incident 
Cassegrain telescope mirror. The orientations of the lines 
represent the orientation of the fast axis of the retardance. _ 
The lengths of the lines signify the magnitude of linear 
retardance, which increases guadratically with pupil 
coordinate. This polarization aberration is called linear 
retardance defocus or quadtard. Figures 28.a and 28.b show the 
effect of guadtard on a uniform lef t-circularly polarized beam 
(a) and on a uniform, vertical-linearly polarized beam (b). The 
dominant effect of the retardance is a coupling of linear into 
circular polarized light and vice versa. These figures are for 
retardances about 40 times greater than the retardance of a 
conventional Cassegrain. 

in general, a Cassegrain telescope displays both these 
polarization aberrations simultaneously, linear polarization 
defocus and linear retardance defocus, with the retardance being 
the larger term. This polarization aberration induces 
polarization coupling which reduces the accuracy of subsequent 
polarimetric measurements. The polarization coupling for a 
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SAMEX-type Cassegrain design but with standard aluminum coatings 
is 3 x 10" 3 , almost two orders of magnitude greater than the 
radiometric accuracy of the optics. Similar amounts of 
polarization accuracy would be associated with relay optics 
utilizing standard antireflection coatings. 

3.3 The SAMEX Design 

It is clear from the discussion of the previous section 
that standard optical designs will not suffice for the SAMEX 
foreoptics because they produce unacceptable levels of 
instrumental polarization. The design of the SAMEX foreoptics 
resulted from the development of a new method for the analysis of 
instrumental polarization based on the theory of "polarization 
aberrations." This theory allows the description of the 
variations of amplitude, phase, and polarization of an optical 
wavefront across the exit pupil of an optical system. Because 
the theory naturally incorporates the polarization properties of 
the thin film coatings on the individual surfaces of the optical 
system, this method integrates the coating design with the lens 
design. In ordinary optical design work, these two phases of the 
design are normally decoupled and pursued separately. 

This unified optical and coating design was performed 
for the SAMEX foreoptics to insure that the SAMEX magnetograph 
will accurately measure the polarization state of incident 
sunlight. To improve the SAMEX polarization accuracy over 
standard designs, various telescope and lens coatings were 
investigated. Most standard reflection-enhancing coatings for 
telescope mirrors were found to be substantially worse than the 


45 


6 

•••••• G: C 

' ■ gc 'cc c " 
'Ggcccgg 

CCCGC 
C C G Co 

o 


1 

3. 

- 4 

C 6 1 

Q > 

1 1 1 

1 

3 - 

i i m 


l ^ 

3 i i 

I 

3 

.031 

1 

0 

G _ 

1 


< 


B. 


Figure 28. Linear retardance effects associated with a 
Cassegrain telescope. The large effect of linear retardance 
a Cassegrain mirror is shown for (a) a beam of uniform left 
circular polarization and for (b) a beam of uniform linear 
(vertical) polarization. The dominant effect of linear 
retardance is a coupling of linear polarization into circular and 
vice versa. 
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bare aluminum. Two coatings were designed which had 
significantly improved polarization performance for a 
telescope. These coatings are described in the section on 
coatings (section 3.6). Figure 29 shows the linear polarization 
(a) and linear retardance (b) associated with a Cassegrain 
telescope when these coatings are utilized and when ordinary 
aluminum coatings are used, as a function of the height of the 
ray entering the telescope, from top to bottom. The amount of 
polarization or retardance for the three coating choices is 
plotted along the y-axis for (1) two aluminum coated mirrors, (2) 
two Q201 coated mirrors, and (3) one Q201 coated mirror and one 
eight-layer enhanced-ref lection coating. The third design (3) ^ 

balances polarization and retardance effects of opposite signs 
between the two mirrors to achieve a polarization performance 

superior to either mirror separately. 

Similar design strategies have been used with the relay 

lenses to achieve significant improvements over conventional 
coatings and to balance the remaining polarization effects. As a 
result, the overall polarization performance for the foreoptics 
has been improved by a factor of ten thousand relative to 
conventional designs. Thus the actual attainment of the 
polarimetric sensitivity determined by the SAMEX requirements 
will presumably be limited only by whatever scattering effects 
are present, and not by the instrumental polarization of the 

optics. 

The method used to achieve this dramatic result - 
polarization aberration theory - is outlined briefly in sections 
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3.4 and 3.5 and expanded upon in Appendix C. A full development 
of the method is given in Chipman (1987). The basic results from 
the method are the second-order polarization aberration 
coefficients which provide a quantitative measure of the 
polarization accuracy of any optical system. Specifically/ these 
coefficients determine the parameter A p , the polarization 

accuracy, defined as the maximum fraction of light which can be 

* 

# * 

coupled into an orthogonal polarization state. It is given in 
terms of the second order polarization aberration coefficients 
P(l, 0,2,2), P(l, 1,1,1), and P(l,2,0,0): 

A p < P( 1 , 2 , 0 , 0 ) 2 + \ P(l,l,l,l) 2 + -j P(l,0,2,2) 2 . 

For the SAMEX magnetograph foreoptics and coating design given 
herein, the value of 

A < 1.4 x 10 7 
P 

was obtained. For standard coatings of aluminum, the 
polarization accuracy is 

A = 2.7 x 1CT 3 
P 

for a Cassegrain telescope alone (no relay lenses). (The second 
order aberration coefficients used in these calculations are 
given in Table 19.) This result for A p means we have achieved 
our design goal: the polarization state of the light from the 


49 


ORIGINAL PAGE IS 
QF. POOR QUALITY 


Sun can be determined without introducing polarizing effects from 
optical elements in the system. We have been able to effectively 
eliminate the problem of induced instrumental polarization m the 

SAMEX foreoptics. 

3,4 Polarization Aberration Theory 


The SAMEX foreoptics are intended to transmit all 

polarization states equally. But all optical interfaces display 

some polarization when used at non-normal incidence. Thus 

polarization is present in all systems at some level. If the 

system is intended to be nonpolarizing, the instrumental 

polarization is often termed 'residual polarization" to signify 

its generally undesirable character. Residual polarization might 

be compared to wavefront aberration because both interfere with 

the measurement of optical fields and reduce the image forming 

potential of the optical system. 

The principal cause of instrumental polarization in the 

optical systems of present solar magnetographs is the 
polarization due to non-normal incidence at the optical 
interfaces'” and coatings. Since each ray takes a different path 
through ^ system with its own angles of incidence and planes of 
incidence, each ray in general experiences a different change in 
its state of polarization. This residual polarization varies 

with wavelength, object coordinates and pupil coordinates. 

, . . _ it ,.:i i hp defined as variations of the 

"Polarization aberrations will b 

. ^ nmi ari ration of an optical wavefront across 

amplitude, phase and polarization 



the exit pupil of an optica! system and the dapendence of these 

variations on wavelength and object coordinate. The polartzatron 

aberrations are extensions of the wavefront aberrations and 

encompass both ampiitude and polarization variations, thus 

providing a more complete characterization of the electrom g 

fields transmitted by an optical system. 

vacuum-deposited thin films are used on most opticai 

surfaces to control the amount of transmission and reflection. 
These thin films are usually less than the wavelength of Ug 
thickness. Seing so very thin, the effect of the films on ray 
paths are accurately modeled by treating the films as having 
parallel surfaces which contour the substrates on which they are 
deposited. Cue to the closely spaced parallel surfaces, t.in 
films have negligible influence on the ray paths through the 
aystem and are generally ignored when simulating a system by ray 
tracing. These coatings principally affect the amplitude and 
polarization of the ray and have much less effect on the optica 
path difference. This division, with the optical surfaces 
governing the ray paths and the thin film coatings governing 
amplitude and transmission, allows the optical system design 
problem to neatly decouple into- two separate problems, lens 
design and coating design. The wavefront performance and image 

• dilated bv a lens designer using a 

quality of the system is calcula Y 

ray tracing optical design program. The amplitude and 
polarization calculations at individual surfaces are performe 

using a thin film design program. 
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This decoupling of optical design and coating design has 
usually worked well. The coatings designed to optimize the 
transmittance or reflectance at an interface have usually reduced 
the amplitude and polarization variations and thus reduced the 
polarization aberrations at the interface as well. For example, 
a quarter-wave magnesium fluoride antireflection coating on glass 
typically reduces reflection losses. at the design wavelength by a 
factor of four, and reduces the instrumental polarization by a 
comparable factor. This fortuitous circumstance has allowed lens 
and coating design to remain relatively uncoupled. Thus 
instrumental polarization was usually ignored as a higher order 
effect. But it is not sufficient to design thin film coatings_in 
isolation from the lens design for the SAMEX magnetograph - the 
demands on amplitude and polarization performance are too great. 

For the designs of the SAMEX system special methods 
have been developed to- calculate the instrumental polarization of 
the SAMEX foreoptics. These methods are described in detail in 
Appendix C. Calculating the instrumental polarization requires 
performing thin film calculations during the ray tracing 
process. This Idea is not new, but Its implementation is complex 
enough to have delayed this obvious integration of these two 
branches of optical design until specifications required it. 

In this new methodology, a Jones polarization matrix is 

calculated for arbitrary optical paths through the optical system 
and includes the effects of all the specified optical coatings on 
the curved optical interfaces. The specific technique used for 
this calculation is the method of -polarization aberrations." 
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Th U is a method analogous to the aberrations of geometrical 
optics (spherical aberration, coma, astigmatism, etc.) except 
th at it encompasses the polarization effects of coatings as well 

as the wavefront aberrations. 

The analysis of the SAMEX optics proceeds In several 

stages and Is summarized In Table 16. first, the optical system 

is designed using the CODE V optical design program to optimize 

the optical design for a spatial resolution with a minimum 

aperture. During this phase the angles of incidence are kep 

small as possible to reduce polarization effects from coatings. 

Second, special thin-film, reflection-enhancing 

designed for the telescope mirrors and special antireflection - 

coatings are designed for the lenses because conventional .coat ng 

designs show significant polarization and retardance 

contributions near normal incidence, these coatings are specie 

designs which reduce the polarization effects at small angles o 

incidence over the wavelength band of the magnetograph. Next, a 

Taylor series is calculated to represent the coating 

performance, finally, the effect of these coatings In the 

magnetograph optical design is calculated using the polarization 

-w a ,-hi« last step produces the parameter A , the 
aberration method; this last seep p 

polarization accuracy. 
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Table 16- S«pe In Instrumental Polarization Simulation 
1 Design optical system using CODE V. 

2 . Design low polarization, high reflectivity telescope coatings 

3. Design low polarization, antiref lective lens coatings. 

4 . Determine Taylor series representation of coating 
performance . 

5. Calculate polarization aberration .coefficients of optical 

system • 

6 . Calculate A p . 

7. Iterate 2, 3, 4, 5 5 6 until satisfactory performance is 

achieved ( A p < ^ )• 


in the following section, the method for deriving the 
polarization aberration coefficients and the polarization 
accuracy a p is outlined in more detail. 

3.S Polarization Aberration Coefficients and Polarization 
Accuracy 

The polarization states of the electromagnetic field 
are described by the complex two-component Jones vector, 3 . 
polarization states and hence the Jones vector are transformed 
when a ray passes through an optical interface which is described 
by the Jones matrix transformation JJ. For the SAME* foreoptics 
we have homogeneous, weakly polarizing optical 

design, for which the transmission coefficients perpendicular, 
t and parallel, t p , are given in terms of a Taylor series 
expansion of the angles of the chief and marginal rays, l c and 
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tsspec t i v© ly • Each individual optical ray path tor the 
rotationally symmetric system can be defined in terms of the 
object height, H, and the pupil coordinants, p and <f> as defined 
in Figure 30. The cascade effect of an optical train is given by 
the product of the individual Jones transformation matrices for 
each element interface. Then the overall polarization of the 
foreoptics is described by the complex Jones transformation 
matrix which is expanded in terms of the ray coordinates 
(p, and H) and the basis matrix set cr ( k ) : 

JJ = E E EE P(k,u,v,w) H U p V cos W ( $ ) o(k) 
k=0 u v w 

where P(k,u,v,w) are the expansion coefficients. The 2 by 2 

matrices, a(k), which define the k th polarization state, are- the 
identity matrix and the Pauli spin matrices. This expansion is 
in the same form as the standard wavefront aberration expansion.- 
These polarization expansion coefficients for an expansion to 
second order (specified by s * 2) in the angles of incidence are 
a function of the total optical transmission, t, the normalized 
secondard order basis set coefficients, d(k,s,q), for each 
individual optical element, q, and the angle of incidence of the 
chief and marginal rays. Hence these coefficients are dependent 
on the characteristics of the optical coatings and the ray 
tracing results. These SAMEX specifications are given in Table 
19a-c . 

The polarization aberration expansion for the radially 
symmetric system of interfaces with isotropic coatings and 
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figure 30. Paraxial coordinate ^^^'em^Th^^axis is the 

jp^cal* ax is^of ^a^ rota t iona lly synme tr ic^opt ical *^ s th|^ u 9^ 

initially travels ™ “Seized by ray coordinates at the _ 

in optical system are char normalized object 

abject and entrance pup • li2ed pu pil radius, and ♦ is the polar 
coordinate, p ls .^ h sure( j counterclockwise from the y ax ^s. 

angle in the P u P^ r ””an ctorllnates in the pupil are x and y. 

The normalized Car ***} a " ^a re also shown. 

The chief and marginal rays are axs* 



nonpolarizing transparent media is given 
angles of incidence by the polarization 


to second order in the 
aberration expansion 


JJ(H,p ,♦) 


t [ a Q + P(l, 2 , 0 , 0 )H 2 a 1 + 

P(l/l/l/l) H P( a 1 cos<> - a 2 sin4>) + 
P(l,0,2,2)p 2 (a 1 cos2«D - <r 2 sin2*)l. 


The amplitude a ljkl of the polarization aherratlon coeftzczent 
P,i (j .K.l> describes the magnitude o£ the polarization aberration 
effects while the phase < ijkl of the coefficient describes the 
magnitude of the retardances. Then, given a specific incident 
polarization state, the polarization state of the exit beam 
be calculated from this matrix. Therefore the amount of 
polarization in the orthogonal state can be determined. This 
orthogonal polarization in the exit beam then determines 

accuracy of the polarization measurement. 

As an example of the coupling of the optics to an 

orthogonal state of polarization, i.e., polarization crosstab, 
consider the following example. The on-axis linear polarization 
and linear retardance of the S^EX foroptics, i.e., the term 
linear defocus, is described by the fourth term in the 
polarization aberration expansion JJ. The instrumental 
polarization function J d (H,p,») for linear defocus is then 


* x[. 0 * PU.0,2,2>» 2 <‘1=°^ ' 

2 (<j^cos2$ - cr 2 sin2$)]p 


0 

= x[a Q + ( a iQ22 


+ 1 ^ 1022 p 
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where t is the amplitude of transmittance of the system down the 
optical axis. The transmittance (t) describes the polarization 
independent reflection and absorption losses associated with the 
ray down the optical axis. P(l,0,2,2) describes the linear 
polarization (a 1022 ) and linear retardance U 1022 ) associated 
with the marginal ray. 

At any given point in the pupil/ the eigenenpolarizations 

1 A 

are linearly polarized light, oriented radially, and 

tangentially J t . This concept is important to the definition of 
polarization accuracy which is to be defined. We shall now 
calculate for circular and linear incident polarize light to 
illustrate of the effect of crosstalk. 

Maximum coupling occurs when the incident light is 
circularly polarized, since circularly polarized light can always 
be decomposed into equal components of J r and J fc everywhere m 
the pupil. The coupling is zero in the center of the pupil 
(where the polarization and retardance vanish) and increases to a 
maximum coupling at p = 1 of 


I c,max (H,1 '* > 3 |Pd»0,2,2) | 2 * a ? 0 22 + 6 1022 

at the edge of the pupil. The net fraction of incident 
circularly polarized light coupled into the orthogonal circularly 
polarized state, I*, is given by the integral over the pupil, 

2 ^ i 2 

jl = Jk f d$ / pdpjp P(l,0,2,2)| , 

c * o 0 

|P(1,0,2,2) | 2 
= J ^ 


I£ we have a Stokes vector o £ W then the amount o£ polarization 
in the orthogonal state -V can be calculated. If -V U the 
maximum amount of crosstalk, we identify this as the polarization 

accuracy. 

For incident linearly or elliptically polarized light. 
Erection of coupled intensity is less because the light is not 
composed of equal fractions of eigenstates. The coupling is 
minimum for incident linearly polarized light, which will be in 
one of the eigenpolarizations along one axis in the pupil and in 
the orthogonal eigenpolarization along the orthogonal axis. Here 
the fraction of coupled energy will be calculated assuming an 
incident polarization state of horizontal linearly polarized 
light A for calculational simplicity; the same fraction is 
coupled for any linearly polarized incident state. The 
orthogonal state of vertical linearly polarized light is 
designated as V. The polarization state transmitted by an 
optical system described by linear polarization defocus for H is 


3(h,p/*) 


tIH + P(1,0,2,2)p 2 (H cos 2 + - V sin 2 *) 1 . 


The fraction of incident A light coupled into V light is equal to 
o 1 * i 2 

t 2 = I r w do / P dp|3(H,p,$) • v| 


- I f 2 ' d* f 1 pdp|P(l.°, 2 , 2)» 2 sin 2 »| 

* 0 o 
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a )P(1>°> 2 > 2 ) 1 f s in 2 2<frd$ / 5 p dp 
* 0 0 

|P(1,0,2,2)| 2 

= j 

Therefore, we defined the polarization accuracy, * p , as the 
maximum fraction of light which can be coupled into orthogonal 
polarization states. The incident polarized state is given by 
the Jones vector, 3. This vector is rotated by the optical 
system and the rotation is given by the Jones matrix, JJ. The 
amount of polarization along the orthogonal state of 
polarization, 3’, of the incident polarization state is given by 
the projection of JJ(J) into <?’• i.e. 

JJ(3) *3 ' . 

The maximum projection is used in defining the polarization 
accuracy. In the above example ij > I c 2 / and only I C 1 is 
retained for the linear defocus term, P(l,0,2,2). 

This value is given by the square of the second order Jones 
matrix and is given in terms of the polarization aberration 
coefficients. Then the polarization accuracy is given in terms 
of the second order polarization abberation coefficients 
P(l,0,2,2) , P(l, 1,1,1)/ and P( 1,2, 0,0): 

A = — J I JJ,(l,p, + ) 1 2 pdpd* 

p IT ' 1 ^ 
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or 


» < |P(1,2,0,0)| 2 *\ |P<1.WM 2 l p ‘ 1 '°' 2 ' 2, | 2 ' 

P " ’ 

where the integrations have been carried out for the squa 
terms and estimated for the crossed terms as being less than or 
equal to the direct product terms. For the SAMEX magnetograph 
design, the value for the polarisation accuracy. 


A < 1.4 x 10 
P 


-7 


is obtained for specially selected optical coatings. These 
second order aberration coefficients are given in Table 19. 
second order coefficients are sufficient since the next order 
that contributes is the fourth order and their polarisation 
effects would be on thejgrder of U p > 2 - I" the following section 
we discuss the process used to design these special optical 

coatings . 
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Table 19-a 


. Table of Total Second Order Polarization 
Aberration Coefficients for the 
SAMEX Magnetograph Given the Individual 
Surface Components of Table 19-b 


Vec tor Quadratic Piston ; 

P(l, 2,0,0) = * d ( 1 , 2 ,q ) i c 2 (q) 

= -7.72x10”^ - T 1.50xl0~ 4 radians 2 


Vector Tilt; 


P(l, 1 , 1 , 1 ) = 2 * d ( 1 , 2 ,q ) i c (q) i m (q) 


= -1.42xl0~ 5 - T 2.70xl0“ 5 radians 4 


Vector Defocus: 


P ( 1 , 0 , 2 , 2 ) = “ d ( 1 , 2 ,q ) i m 2 (q) 


s -1.91xl0“ 5 + T 1.80xl0" 4 radians 2 
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Contributions 


Table 19-b. 


Surface by Surface 
Given Surfaces as 


Polarization Aberration 
Defined in Table 19-c 




Vector 


Del 


d( 1 » ^ 


ocus 

■(q) 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


8 . 3e-10 
-l.le-8 
-4 . 7e-7 
1 . le-5 
-4 . 9e-6 
-1 . 3e-6 
-5 . le-6 
-2 . 3e-5 
-1 . 3e-5 
-l.le-8 
-5 . 7e-7 
-l.le-6 


3 . 4e-7 
-1 .3e-6 
-1 . le-6 
-2 . 6e-5 
-1 . 2e-5 
0 

-1 . 2e-5 
-6 . 6e-5 
-3 . le-5 
l.le-6 
0 

-2.7e-6 


-3 . 9e-8 
3 . 9e-7 
4 . 4e-7 
-3 . le-6 
3 . 6e-7 
-2 . 3e-6 
-1 . 4e-6 
-2 . 3e-6 
-1 . 5e-6 
5 . le-8 
1 . 2e-6 
1 . 3e-6 


-3.9e-5 
4 . 6e-5 
5.8e-7 
-7 . 5e-6 
8 . 5e-7 
0 

-3 . 4e-6 
-5 . 5e-6 
-3 . 6e-6 
-4 .9e-6 
0 

3. le-6 


1 .3e-6 
-1 . 4e-5 
-1 . 3e-7 
-9 . le-7 
-2. 6e-8 
-3 . 8e-7 
-4. le-7 
-1 . 9e-7 
-1 . 7e-7 
-2 . 3e-7 
-2.7e-6 
-1 . 5e-6 


1 . 3e-3 
-1 . 7e-3 
-3 . Oe-7 
-2 . 2e-6 
-6. le-8 
0 

-9.3e-7 
-4 . 6e-7 
-4 . 2e-7 
-2 . 2e-5 

0 

-3 . 6e-6 


Note: The surface numbers S(i) are 
q by S ( i ) = q + 2, i.e. q-1 

prefilter is not considered. 


related to the surface 
is the primary mirror and J h f% xl0 -10 
The notation 8.3e-10 means 8.3x10 
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Table 19-c • 


The Surface Definitions for the Polarization Aberration 
Coefficient Calculations 


Surface Optical 
q Coating 


Paraxial Angle 


Second Order 
Polarization 


** 


Re 


Im 


Magnitude of 
2nd Order 
Polarization 


1 

2 

3 

4 

5 

6 
7 
3 
9 

10 

11 

12 


refl201 
layers 8onAl 
v!0c5256 


kf9fk5 

vl0c5256 


vl0c525 

kf9fk5 

vcl0c5256 


0.00134 

0.00236 

0.03862 

0.18664 

0.12576 

-0.30317 

-0.12743 

-0.29808 

-0.20326 

-0.00703 

-0.06243 

-0.5980 


-0.06250 

-0.08519 

-0.02031 

0.05388 

-0.00908 

-0.05091 

-0.03622 

-0.02495 

-0.2363 

0.03169 

0.13588 

0.06974 


2.82e-7* 
-1 . 20e-6 
-1.90e-7 


-8.91e-8 

-1.90e-7 


-1.38e-7 
-8. Ole— 8 
-1.90e-7 


2.84e-4 
-1.39e-4 
-4 . 53e-7 


0.0 

-4.53e-7 


1 . 33e-5 

0.0 

-4.53e-7 


2.8e-4 

1.4e-4 

4.9e-7 


8.9e-8 
4 . 9e-7 


1.3e-5 
8.9e-8 
4.9e— 7 


*(2.82e-7 = 2.82xl0“ 7 ) 

**The rate of .linear polarization (Re) and retardance (Im) per angle of incidence 
squared (deg -2 ) for the specified optical coating. 
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3.6 Coating Designs 

A design study was performed to find coatings especially 
suitable for use in the SAMBX magnetograph. First, many 
conventional coatings were examined and found to have 
polarization and retardance properties that make them 
unacceptable for use in the SAMEX magnetograph. Alternative 
coating designs were devised with improved polarization 
performance and the best of these were specified for the 
instrumental polarization simulation. Coatings were designe 
optimized and analyzed using two thin film coating programs: FTG 
Software's Filmstar and Prof. A. Macleod's (University of 
Arizona) FILMS. 

The requirements for the SAMEX magnetograph coatings ace the 
following : 

1) They must have extremely low polarization properties 
near normal incidence so that they introduce a minimum of 
polarization or retardation into the optical system. 

Specifically, the differences between the s and p amplitude 
coefficients for reflection or transmission for less than 10* 
angle of incidence should be 0.01% and the phase difference 

should be less than 0.05°. 

2 , The system's transmission must not be compromised by the 
coatings: reflective coatings should reflect >09.9% and 

transmission coatings transmit > 99.7%. 

f,nH,rahie and not be excessively 
3) They must be manuf acturaDie 


costly . 
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4) The amount of scattered light should be small. 

5) They should have reasonable tolerances for fabrication. 

Several common reflecting coatings for telescopes and 
antireflecting coatings for lenses were analyzed and found to be 
inadequate for the SAMEX magnetograph. For example, telescopes 
are frequently coated with bare aluminum coatings (Al) or 
quarter-wave stacks on top of aluminum. These coatings have very 
little linear polarization near normal incidence but substantia 
amounts of linear retardance. The common antiref lection coatings 
for lenses (quarter-wave MgF 2 and broadband antiref lective 
coatings) have very little retardance near normal incidence but _ 

substantial linear polarization. 

Three coatings (see Table 20) have been selected from our 
coating design study for inclusion in the SAMEX magnetograph 
design. The polarization effect of these coatings has been used 
in the instrumental polarization simulation. These coatings 
demonstrate that coatings with the necessary low polarization 
performance can be designed. The present coatings are somewhat 
sensitive to fabrication errors and are very wavelength 
sensitive. The designs would benefit from further work to 
understand why they work, knowledge which should lead to fully 
optimized designs with improved manufacturing tolerances. 

Table 20 contains the description of the three coatings 
specified. In addition, the optical system contains two cemented 
interfaces, between kf9 glass and fikS glass, specified kf9fk5. 

The polarization of this interface, although negligible, was 
included in the instrumental polarization calculation. 
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p fdiire 31. "Optical properties of the interference coating Q201. 
•he phase and reflectance properties of the interference coating 
)201 P designed for the SAMEX telescope mirrors are shown as func- 
:ions of the angle of incidence (A, on the left) and wavelength 
on the right). ( A 1 ) The phase of a reflected beam is shown 
as'a function of the angle of incidence. The s and p components 
3 f an incident wave of wavelength 5250 A are shown as one curve. 
(A2) The difference in phase of the s and p components shown in 
/ Ai ) is clotted as a function of the angle of incidence. (A3) 
rhe % reflectance of a 5250 A reflected wave is shown as a func- 
tion of the angle of incidence. (Bl) The phase of the reflected 
wave is plotted as a function of wavelength. ( B2) 
reflectance of a wave incident at an angle of 2 is shown as a 
function of wavelength. (B3) The percent reflectance of a wave 
incident at an angle of 6« is shown as a function of wavelength. 
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The coating 0201 Is a 201-layer, quarter-wave enhanced 
reflection coating specified for the primary mirror. The 
philosophy behind this design is that the polarization effects 
arise because of large refractive index differences between 
layers. This coating uses two materials with a small index 
difference, fused silica at n-1.45 and evaporated glass at 
n =l . 52 . To achieve the desired reflection performance with 
smaller index differences requires more layers, 201 in this 
ease. The materials chosen are both amorphous and thus quite 
suitable for coatings with large numbers of layers, since they 
don’t display the microcrystalline growth patterns which lead to 
the unacceptable scattering and inhomgeneity associated with _ 
crystalline materials. Figure 31 shows performance curves for 
Q201 both as a function of angle of incidence and wavelength. 
Figure 31..-1 shows the absolute phase change on reflection for 
the s and p components. The quadratic portion of these curves is 
the defocus introduced by the coating. The parameter of greatest 
interest is the retardance, the difference between the s and p 
phase changes; it is plotted separately in 31.a-2. This 
retardance is below 0.1* over the angles of incidence of 
importance. Figure 31.a-3 shows the s and p intensity 
reflectance as a function of angle of incidence. The difference 
between the s and p reflectance is the linear polarization, here 

less than 0.01%. 

The coating specified for the secondary, LSonAl, is a more 
conventional 8-layer, reflection-enhancing coating specially 
optimized to complement the Q201 coating. This coating has the 
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opposite sign on both its linear polarization and retardance 
relative to 0201. Thus the small residual polarizations of these 
two coatings tend to cancel, resulting in a nearly polarization 
free Cassegrain design. 

The coatings specified for the lens surfaces are a two-layer 
coating of the V coating family of designs. A region of solution 


was found with the remarkable property that the linear 
polarization and linear retardance both changed signs within 50 
nm in wavelength of each other. Thus, by varying the thicknesses 
of the layers, low polarization coatings with several different 
useful properties were obtained. The two scalings used in the 
design were 1.0 for the coating vl0c525 and 1.06 for vcl05256. 
These coatings have by far the lowest polarization effects of any 
lens coatings investigated. Fabricating and testing samples of 
these coatings should be. conducted in the next study phase. 


Table 20. Coatings Specified for the Magnetograph. 


Name 


Number of Polarization Reflectance or 
Layers Magnitude Transmission 


Coating 

Design 


V10C525 


Q201 

.82L) 


L8onAl 

,25L)2 



Reflective 


1.4e-4 


99.7% 

( .25H .25L) 100 .254H 

H=l. 52 L=1.45 

2.3e-4 

99.98% ( .062H 


H=2. 15 L=1 . 38 

99.65% 

. 38H . 29L . 26H 

• 25L ( .25H 


H=2. 38 L=1 .38 
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have 


A11 the polarization analyses of thin films performed here 
assumed ideal, thin-film structures, materials which are 
uniform, isotropic, homogeneous and free from scattering. This 
is a sufficient assumption for the design and analysis of the 
thin films. Actual thin films possess extremely complex 
microstructure which contributes to scattering, inhomogeneity and 
anisotropy. The typical thin film microstructure is an array 
columns growing up out of the substrate. This causes nonideal 
coating properties which will impair the performance of the 
system at some level* 

The performance of the SAMEX magnetograph will be enhanced, 
by using the finest coatings available, which will probably be 
deposited using ion-assisted thin film deposition. By bombarding 
the growing thin film with energetic ions, usually helium or 
argon, the growing thin film structure is pressed down and 
becomes denser (Martin, 1986). This disrupts the growth of the 
lat ,e columnar structures which are responsible for much of the 
scattering and anisotropy. The resulting films have lower 
scatter and greater uniformity. The ion-assisted thin film 
deposition technique is still new and not yet widespread in its 

use . 

As part of the preliminary efforts leading to the actual 
fabrication of the SAMEX magnetograph, test coatings need to be 

prepared and coating vendors qualified. Several coating 

, , , -Anhrarted to produce small (2x2 inch) 
manufacturers should be contracted P 

. , coat inqs should b© 

samples of the specified coatings. These coat g 
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tested for spectral transmission, polarization performance and 
scatter (bidirectional reflection measurements) as part of the 
final design and vendor qualification. The Marshall Space Flight 
Center and the University of Alabama in Huntsville Center for 
Applied Optics have the expertise for such coating 

characterization. These data would then be used to calculate the 
impact of the stray light coming from the coating scatter on the 
magnetic field accuracy of the magnetograph. 

3.7 Further Sources of Instrumental Polarization 

Two effects that could degrade the polarization performance 
of the SAMEX optics are coating anistropy and scattered light.— 
Coating anisotropy is the variation of refractive index with 
direction in the coating. Coating anisotropy can be measured 
ellipsometrically. It frequently occurs where coatings have' been 
deposited at non-normal incidence. Most coatings have a columnar 
microstructure which usually grows out of the substrate toward 
the source. For, a nonnormal deposition angle, the resulting 
coatings have builtin birefringence. For the magnetograph 
optical design, this additional coating-induced retardance is a 
problem which must be held to acceptable levels in the final 
coatings. A level of anisotropy below a root sum squared 
birefringence of 0.0001 radians rms per coating is highly 
desirable . 

Scattering depolarizes light: the scattered light is random 
and carries less information about its original polarization 
state. Coatings, because of their detailed microstructure, may 
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display substantial scattering. Considerable effort is being 
devoted on many approaches to producing optics with reduced 
scatter, including the ion-assisted deposition coatings already 
mentioned. Scattering in the coatings of the foreoptics will 
probably be the limiting factor in the accuracy of the 
magnetograph, now that the instrumental polarisation of the 
Cassegrain telescope has been reduced by orders of magnitude. 


3.8 Optical Tolerances 

An optical tolerance study of the magnetograph is outside 
the scope of this study. However we would Ilk. to point out the 
areas of special study that are needed for a flight instrument,. 
The individual optical components need to be qualified for 


spaceflight. This would include operation under a vacuum for an 
extended period of time. The effects of particle radiation over 
the lifetime of the spacecraft must be resolved. The mechanical 
robustness of the optical system is to be defined. This would 
include the effect of thermal drift on the optical alignment. 


The general effects of mechanical fitting and sensitivity due to 
fabrication errors need to be addressed, including surface 
figure, tilts, decentrations, and optical coating variations. 

The overall alignment and calibration procedures are to be 
specified before the optical design is completed (Yoder, 1986). 
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4 . CONCLUDING REMARKS 


In this study we conclude theoretically, that the electromagnetic 
soliton in a current sheet could trigger the resistive instability for an 
uniform magnetic field configuration, which eventually turns in an 
eruptive instability at the onset of the magnetic field reconnection. 
Therefore, we may consider this physical scenario which could be a 
candidate for solar flare triggering mechanism. Experimentally, we 
presented an analysis of instrumental polarization of imaging optics. 
This optical system cold accurately measure the polarization state of 
sunlight in a narrow spectral bandwidth over the field view of an active 
region to make an accurate determination of the magnetic field in the 
region. 
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appendix 


POLARIZATION ABERRATION THEORY 


Introduction 

This appendix describe the method which have been developed 
to calculate- the instrumental polarization ot the SAMEX optics. 

A full development of the method is contained in "Polarization 
Aberrations" by Chipman (1987). The calculation requires 
performing thin film calculations during the optical design 
process to determine the two parts (eigenvectors) of an optical 
beam, as a function of the object and pupil coodinates. 

The mathematical method of treating instrumental 
polarization will be first discussed in terms of the Jones and C- 
vectors and then the polarization aberration expansion will be 

derived. 


The Jones Matrix and C Vector for the Characterization of 
Polarization 

The most efficient mathematical method for treating the 
SAMEX instrumental polarization is the Jones calculus. The 
Mueller calculus is a more difficult representation which 
includes optical depolarization (scatters, properties. Such weaX 
depolarization effects are more readily handled by experiential 
measurements. The Jones calculus (Jones 1941a-c , !942 , 1947 , Azzam 
and Bashara 1977, Theocaris and Gdoutos) is a mathematical 
formalism for treating problems involving the description of 
polarized light and polarizers which uses the Jones vector for 


the description of polarized light and the Jones matrix to 
characterize the polarizing properties of an optical element. 
The details of the following discussion are given by Chipman 
(1987). The elements of the Jones matrix and the C vector 
characterization are outlined first. 

Definition of the Jones Vector in terms of the Electric Field 
Amplitudes 

The Jones vector expression for the a quasi-monochromatic 
plane wave propagating parallel to the z axis with electric fie 

amplitude 

E { t ) = E x (t) + E y (t) , 


where , 


E (t) = n x E 0f . 


(t) cost(kz-wt) + e x l ' 


and , 


y * n y E 0,y 


(t) cos [{kz-wt) + e y l » 


where S* and S y are direction unit vectors In the * and y 
direction with the light propagation in the z direction. There 
are 4 parameters E 0fX . E 0 , y , e x , and e y defining the wave beside 

the wavelength. 

The time dependent Jones vector is defined in terms of the 


electric field amplitudes as ; 


The components of J(t) are the instantaneous components of E(t). 

The normalized Jones vector J is a time independent 
normalized vector where all the vector components of J(t) have 
been divided by the incident electric field amplitude, 

J = J(t)/E 0 (t) . 

The normalized Jones vector is referred to as "the Jones vector"^ 
unless otherwise stated. Knowledge of J and E Q provides all the 
information necessary to reconstruct E(t) to within a constant 
phase factor. 

Table C-l lists the Jones- vectors for the most common 
polarization states: horizontal linear, vertical linear, +45 -- 
degrees linear, -45 degrees linear, rightf circular and left 
circular polarized light . These vectors can be multiplied by an 
arbitrary phase factor without changing the polarization ellipse 
of the light; it only changes the absolute phase. 

Having established the vector which defines the polarization 
state we now consider the matrices which represent the 
polarization effect of the optical elements and allow 
polarization calculations to be performed. 
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Table C-l. The Basic Jones Vector Representation for Linear 
and Circular Polarized Light 


r inaAT" PolaiTiz6(i Liaht 



Horizontal Vertical +45 Degrees 

S| l!l 3 111 

C i rcu lar Polarized Light 


-45 Degrees 

s l-ll 


Right Circular 

• i i 

si 1 I 


-1 


Left Circular 
s I 1 
i 


where s = 


✓ 2 


in his original paper, Jones (1941a) shows that the 
relationship between the Jones vector incident on a polarizer, J, 
and that transmitted or reflected by a polarization element, J', 
can always be related by a matrix, the Jones matrix, JJ. Only 
certain transformations of the field components are allowed, 
those describable by a matrix. Thus the fundamental relationship 
between the vector components of the electromagnetic fields 
before and after a polarizing element is, 

J* a JJ J 


, . tt is a two by two matrix with complex 

The Jones matrix, JJ , is a uwu 

elements , 


JJ 


’ j ( 1,1) 
3 ( 2 , 1 ) 


j(l,2) 

j<2,2) 



where j(k,l) * a(k,l) + i b(k,l) . 

Thus the Jones matrix has eight degrees of freedom. Thus 
there are eight different forms of polarization behavior. These 
eight forms are listed in Table 2E. Every Jones matrix 
corresponds to a physically realizable polarizer. 

Sequences of Polarizers 

The Jones matrix associated with an optical ray path through 
a sequence of polarization elements is just the matrix product of 
the Jones matrices for the individual polarizers. If an optical 
ray traverses a series of elements, 1, 2, ... Q, and the Jones 
matrices appropriate to that ray for each element are, JJ(l)r 
j j ( 2 ) , ... J J ( Q ) , then the Jones matrix describing the 

polarization properties of the system along this ray path is 
given by the matrix product, 

JJ = JJ(Q) ... JJ ( 2 ) JJ(1) 

Since the Jones matrix of an optical element is dependent 
upon the wavelength, angle of incidence, orientation, and path 
through the element, each ray in each wavelength will usually 
have a different Jones matrix. Only if a collimated 
monochromatic beam through a series of planar optical interfaces 
can be assumed, then a single Jones matrix can be written for the 
entire cross section of the beam. This is the case in the 
polarimeter section of the magnetograph but not in the 
foreoptics . 


Coordinate System 

The Jones matrix is defined relative to an arbitrary x and y 
coordinate system. Since these coordinates have been defined for 
the Jones vector, the coordinate system of the Jones matrix is 
defined in terms of the Jones vector coordinates. 

It is often desirable to align the Jones vectors coordinates 
with the s and p planes of an optical interface. Only for plane 
surfaces does the orientation of the s and p planes remain fixed 
across the surface. For nonplanar surfaces, it is necessary to 
maintain two sets of coordinates, the global x and y coordinates 
with respect to which the Jones matrix is defined, and a local s_ 
and p coordinate for each individual point on the interface. The 
local s and p coordinate system will have its x’ and y 1 axes 
aligned with the local s and p planes of the surface for the 
evaluation of Jones matrices at given surface coordinates. Then, 
these local Jones matrices -ill be rotated to bring all the 
matrices into the global coordinate system. 

Pauli Spin Matrix Basis and the C Vector 

The Pauli spin matrices form a most useful basis for 
interpreting the Jones matrix -space- and define a basis set for 
the JJ matrix. The identity matrix, «(0> and the Pauli spin 
matrices, »(1), v(2), and <J(3>, are defined by: 



<y( 0 ) 


<j(1) 

o(2) 


o(3) 

1 

0 


1 0_1 


0 1 


0 

-i 

0 

1 


0 -1 


1 0 


i 

0 


An arbitrary Jones matrix will be expressed as, 

3 

jj » Z c ( k ) o ( k ) . 
k = 0 

The c's are formed into a four element complex vector, the 
"C vector". The C vector expression 


C - [ c(0), c(l), c(2), c(3) ] , 


an equivalent representation of the Jones matrix, 

JJ = c(0) a ( 0 ) + c(l) o(2) + c(2) o(2) + c(3) o(3) 


c ( 0 ) + c( 1) 
c ( 2 ) + ic(3) 


c( 2) - ic ( 3 ) 
c( 0) - c(l) 


When necessary, » end ♦ refer to the amplitude and phase portions 
of the C vector elements, 


C 


[ p ( 0 ) exp ( 7<fr ( 0 ) ) # p( 1 )exp( i ♦( 1 ) ) 


p(2)exp(x4>(2) ) , p ( 3 ) ) exp ( Tp ( 3 ) ) ] 


The elements of C are related to 
equations : 


the Jones matrix elements by the 


c ( 0 ) » (j(l#l) + j(2,2))/2 , c ( 1 ) = ( j (1 ' 1) 

c(l) = + j ( 2 , 1 ) )/2 , c(3) = (j(W2) 

The elements of the Jones matrix are related to the elements of C 
by the equations: 

j ( l,l) = c ( 0 ) + c(l) , j ( 1 ' 2 ) = c(2) - i c(l) ' 

j ( 2, 1) - c( 2) + T c ( 3 ) , 3(2,2) = c(0) - c(l) • 

The C vector, like the Jones matrix, has eight degrees of 
freedom. Table C-2 contains a description of the meaning of th- 
real and imaginary parts of the c vector elements. 


- j(2,2))/2 , 

- j ( 2 , 1 ) )/( -2T) . 


TABLE C-2. 


Interpretation of the C Vector Elements 



Matrix Coefficient 


Meaning 


o ( 0 ) P(0) 

o( 0) *(0) 


Amplitude 

Phase 


Absorption 

Phase 


a ( 1 ) p(D 

a ( 1 ) ♦(!) 


Amplitude 

Phase 


Linear Polarization along Axes 
Linear Retardance along Axes 


<j(2) p ( 2 ) 

o ( 2 ) 4> ( 2 ) 


Amplitude 

Phase 


Linear Polarization, 45 deg 
Linear Retardance , 45 deg 


o(3) 

o(3) 


p ( 3 ) 
♦ (3) 


Amplitude 

Phase 


Circular Polarization 
Circular Retardance 


Matrix and C Vectors 


The Jones Matrix and C Vectors for Specific Polarizers 

Tables of Jones matrices for various polarizers are found in 
A«am and Bashara ( 1977 , Section 2 . 2 . 3 ), Hecht and Zajac ( 1974 , 
Table 8.6), Shurclif f < 1962 , Appendix 2 ), and Theocaris and 
Gdoutos ( 1979 , Table 4 . 1 ). A table listing of the Jones matrices 
and C vectors for the most common polarizers and retarders is 
given in Chipman ( 1987 , Table 6 ). 


The Meaning of the Coefficients of the C Vector 

The primary reason for the introduction of the C vector is 

to simplify the representation of polarizers. Bach of the 
elements of the C vector represents a specific type of polarizer- 

behavior. 

The real parts of the C vector all correspond to amplitude 
effects, absorption and dichroism. The phase portion of the C 
vector represent phase effects, propagation and birefringence. 

The first element, c( 0 )-» ( 0 >exp( i 4 ( 0 ) ) , is the coefficient of the 
identity matrix. Thus it must represent effects that are 
polarization state independent; these are amplitude and phase. 

The last element, c ( 3 )-»( 3>exp( i 4< 3 ) ) , multiplies the spin matrix 
. t3 , which is rotation invariant. Thus the c(3, term represents 
the circular polarization effects, .(31 describes circular 
polarization or circular dichroism and 4(3, describes circular 
retardance or circular birefringence. The remaining two 
elements, c(l> and c(2), represent linear polarization and linear 
retardance. Linear terms require two degrees of freedom: 
magnitude and orientation. Thus, pU> and .(2, characterize 


linear polarization or linear dichroism, p(l) in the 0 degrees 
and 90 degrees directions, p(2) in the + and -45 degrees 
directions. Likewise, $(1) and *(2) characterize linear 
retardance or linear birefringence. 

Rotated Polarizers 

If a polarizer with Jones matrix JJ is rotated through an 
angle 9 (positive if counterclockwise), the Jones matrix becomes 

JJ' (9) = RO) JJ R(-«)« 


The R( 9 ) ' s are the Jones rotation matrices: 

cos ( 9 ) -sin(9) 

R(9) = 

sin ( 9 ) cos ( 9 ) 

The Jones rotation matrices obey the relations, 
R(a) R(b) * R(b) R(a) * R(a+b) 

and, 

R(aJ R( -a.) a a(0) . 

The identity matrix is invariant under rotation; 
R ( 9 ) a ( 0 ) R(-9) = o(0). 


Under rotation, o( 1) and a( 2) couple into each other 


r( 9 ) <j(l) R(-9) = a( 1 )cos ( 29 ) + <x(2)sin(29) 
r ( 9 ) cr ( 2 ) R(-9) = -o( 2) sin( 29) + a( 2 ) cos ( 29 ) . 

a(3) is invariant under rotation; R(0) a(3) R(-9) = a(3). 

Having established the polarization calculus which describes the 
polarization optics, we now apply these matrices to the SAMEX 
optical system. 

Instrumental Polarization 

Two types of polarization calculations can be performed for 
the SAMEX magnetograph: instrumental polarization and transmitted 

light polarization. The first is the calculation of the 
polarization associated with ray paths through an optical 
system. This determines the instrumental polarization function 
as a function of pupil coordinates for a specified object 
point. The other type of calculation determines the state of 
polarization, such as a Jones vector, transmitted by the system 
along a given ray path for a specified input polarization 
state. By iterating this process, the Jones vector as a function 
of position in the exit pupil is calculated. This Appendix deals 
only with the instrumental polarization calculation. Once the 
instrumental polarization function for the system is known, the 
transmitted Jones vectors are readily determined for all input 
polarization states. 


Polarizers are optical elements which divide an optical 
into two parts (Jones vector) and transmit those parts with a 
different transmission coefficient and a different phase, 
two parts of the beam are referred to the eigenvectors or by the 
more descriptive term, "eigenpolarizations . " The two 
eigenpolar izat ions are orthogonally polarized and are transmitted 
by the polarizer with no alteration of their polarization states; 

only the intensity and phase changes. 

The word polarizer will be used to refer to both polarizers, 

such as dichroic or prism type, which have a different 
transmittance for the two eigenpolarizations, and retarders, 
which have equal transmittance but a different phase change for 
the polarizations. -Polarized Light- by Shurcliff ( 1962) is the 
standard reference on the types of polarizers, their definitions, 
parameters and properties. 


Transparent Systems 

The SAMEX foreoptics are highly transparent and weakly 
polarizing and the following calculational method is optimized 
for this case. The ideal Jones matrix for a ray through a 
transparent nonpolarizing system is 


JJ( ideal ) 


exp (Td) 



where d is the optical path length for the ray in radians. (The 
complex value (-IJ^i. denote by ~> • The Jones matrix operation 
on the Jones vector, which is composed of the two orthogonal 



amplitude components, define the transmitted state of the ray. 
The above ideal Jones matrix is the identity matrix, which 
signifies that the system has no absorption or polarisation. 
Since this is the desired form of the Jones matrix for the SAMEX 
foreoptics, the approach developed here obtains the instrumental 
polarisation function as a Taylor series of the system Jones 
matrix in the ray coordinates about JJ( ideal). This approach is 
easily modified for systems which are not highly transparent or 
which contain strong polarisers by performing the Taylor series 
about the Jones matrix for the ray down the optical axis. This 
work deals primarily with this simpler version of the problem, 
transparent systems, to streamline the SAMEX calculations. 


S-P Coordinates 

To handle problems involving light at nonnormal incidence at 
curved optical surfaces, it is necessary to maintain two separate 
coordinate systems: x-y coordinates and s-p coordinates. The x-y 
coordinates are the global x, y and z coordinate system used to 
describe the optical system. The optical axis of radially 
symmetric optical systems coincides with the z axis. 

The s-p coordinates are used to perform polarization 
calculations with the SAMEX coatings which are angle of incidence 

dependent . Most freguently, the functional form of the interface 

mmo^inates. Thus , the rational for 

-polarization is known in s-p 

►Kat- i-voicallV/ the Jones matrix for a 
using s— p coordinates is that/ yp 

. . i i ni-arface will be calculated in the s-p 
ray at an optical interrace 

Then it will be rotated into the system x-y 


coordinates . 
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coordinates. Once all the Jones matrices for the ray at all 
surfaces have been rotated into x-y coordinates, they can be 
cascaded to give the instrumental polarization along that ray 

path in the system x-y coordinates. 

The s-p coordinates are based on the concept of the s and p 
planes. Consider light with unit wave vector k (normalized to 
one) incident at a surface with normal n. The plane of 
incidence, or "p-plane" is the plane which contains k and n. The 
plane perpendicular to the plane of incidence which contains k is 
the "s-plane" . Two unit vectors perpendicular to k are defined 
to form an orthonormal basis for this coordinate system. 

Instrumental Polarization 

All optical elements display variation of their Jones 
matrices as the angle of incidence changes. Further, this always 
involves more than just a variation in the intensity and phase of 
the light; it also involves polarization and retardance. A fine 
optical element used in a transparent system will not display 
polarization effects at normal incidence; it may show some 
absorption, reflection loss or phase shift, but not polarization 

or retardance. 

The Jones matrix can be decomposed and can be expressed as 


JJ 


3 

E c(k) a(k) , 
k=0 


1 0 


1 0 


C o I 1 + 01 ( o -l 5 + 


0 1 0-1 

: 2 ( ) ■+ c 3 L J 

2 1 0 t o 



where a(k) are the identity matrix and the Pauli spin matrices 
which describe the specific polarization state. The c(k) defined 
a vector which then described the polarization properties of the 
element. At an angle of incidence i, the C vector will have the 

form 


C| = (c(0,i),c(l,i),c(2,i),c(3,i)) = 


(p(0 


fi , e t»(°,i), „(2, tie 7 *' 2 ' 1 ', 


where each component has an amplitude (p) and a phase ( ♦ ) part. 
The functional dependences of the C vector coefficients are 
calculated from the Fresnel equations and coating equations for 
the interface. 


Weak Polarizers 

A weak polarization element is defined as a polarizer having 
a C vector such that for some range of i: 


| c ( Q , i ) | 2 >> ( | c ( 1 , i ) | 2 + | c ( 2 , i ) | 2 + | c ( 3 , i ) | 2 ) 


v 2 


Such a polarization element transmits light in a polarization 
state similar to the incident state with only weak coupling into 
other polarization states. The polarization behavior is 
dominated by transmission with only traces of polarization or 
retardance. Any polarization present is at the few percent level 
or less, such that any linearly polarized incident beam has a 
transmission coefficient which varies a few percent or less with 



orientation . 


Similarly, the retardation is degrees or less, far 
less than a quarter wave (*/2) retarder. Near normal incidence/ 
metals in reflection (e.g. telescope mirrors) and dielectric 
refracting interfaces (e.g. relay lens) are weak polarizers. In 
addition, near normal incidence, anti-reflection coated lenses 
used in transmission and metals with reflection enhancing 
coatings are typically weak polarizers for wavelengths near the 
thin film design wavelength. 


Amplitude Transmission Relations 

The amplitude transmission equations for an interface are 
the equations which relate the amplitude and phase of the 
electric fields, E, at an interface. The most general amplitude 
transmission equations for a nonscattering linear interface are: 


“ E' (S) ' 
E* (p) <fc 


'a(ss) E ( s ) + a(ps) E(p) 
a ( sp ) E(s) + a(pp ) E(p) 



where for this section, the plane of incidence will be aligned 
with the y axis. This equation is equivalent to the Jones matrix 

equation. 


" E* (S) 1 


a(ss ) 

a ( sp M 

r 

E(s) 

E' (p) 

1 = 

a (ps ) 

a(pp) J 

E(p) 


For interfaces whose eigenpolarizations are linear polarized 
light oriented parallel and perpendicular to the plane of 


\1 


incidence, the transfer of energy across the interface is 
separable into two uncoupled components which can be written in 

the form: 


E'(s) = a ( s ) E ( s ) = p(s) exp( i<)>(s) ) E(s) 

E’(p) a a(p) E(p) = P(p) exp(iMp)) E(p). 

The amplitude transmission coefficients a(s) and a(p), or in 
polar coordinates, p(s), Ms), P(p), and ♦(?), are determined by 
the Fresnel equations for the interface. The type of energy 
transfer equation, where the s and p equations are separable, is 
a "separable amplitude transmission relation." Only polarization^ 
elements with linearly polarized light as the eigenpolarizat ions 
have the energy transfer equations in. the separable form. This 
includes all the elements and coatings used in the SAMEX 
foreoptics. 

The separable amplitude transmission relations correspond to 
a diagonal Jones matrix in s-p coordinates. The Jones matrix and 
C vector for an amplitude transmission interface in s-p 
coordinates are: 


JJ(i) = 


f a (s , i ) 0 "1 


0 a(p, i ) 


and , 
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C * 1/2 [ a(s, i)+a(p, i) , a( s , i )-a (p, i ) , 0, 0 ]. 

Taylor Series Representation of SAMEX Coatings 

In optical aberration theory, expressions for the optical 
path length of ray segments through the the optical system are 
obtained by performing a Taylor series expansion on Snells law, 
the law of reflection and the grating equation, to obtain 
expressions for the optical path length as a power series 
expansion in the ray coordinates. Thus Snells law, 

n sin i 3 n' sin i' , 
is rewritten for i' as, 

i' = arcsin[ (n/n') sin i ] 

k 

or 

i* = (n/n*) i + [(n/n') 3 ) - (n/n')U 3 /5 + 0(i 5 ). 

The polarization aberrations are generated in an analogous 
fashion. To obtain the variation of the Jones matrix in the exit 
pupil of a system, the appropriate coating equations are required 
in Taylor series form. For radially symmetric optical systems, 
expansions in the angle of incidence about normal incidence are 

used . 

An isotropic interface appears unchanged as it is rotated 
about the surface normal. Ideally, and to a first approximation, 
the SAMEX coatings are isotropic. So, for any isotropic 



interfaces the Fresnel equations are even functions since the 
surface, does not distinguish between angles of incidence of >1 

and -i • 

An even function contains only even terms in its Taylor 
series expansion about the origin. Thus, the Taylor series 
representations of the coating equations has the form, where f(l) 
is a reflection or transmission coefficient, of 


f ( i ) = €q * f 2 i2 + f 4 


i^ + . • • f 


where 


_ . . d n f ( i) 

f n - ^ n! Tpr- 


i=0 


For weakly polarizing interfaces described by amplitude 
transmittance relations, the Taylor series forms of the Jones 
matrix and C vector are calculated as follows. First, the Taylor 
series expansion is determined for the amplitude transmission 

relations : 


( S ) i 2 + a 4 (s) i + • • • / 


,•4 


a(s, i) 3 a 0 (s) + a 2 

a(p, 1) - a 0 (p) ♦ a 2 (p) i 2 + »4 < P> i ' + 

ahou t i=0 in s-p coordinates 
Then, the Taylor series expansion a 

for the Jones matrix is, 



JJ(i) 


a 0 ( s ) 0 

0 a 0 (p) 


+ i 2 


a 2 ( s ) 0 

0 a 2 (p) 


, corresponding C vector expansion in s-p coordinates is 
c = [ c(0,0)+c(0,2)i 2 + ... /■ c(l#0)+c(l»2)i + ••• ' 01 


where 


the nth order c-vector component is given by 


c(0,n) = 1/2 (a n (s) + a n (s)) 

c ( 1 » n) = 1/2 (a n (s) - a n (s)). 

X matrix equation to calculate c(O.n) and ed.n, from thin film " 
program results is given by Chipman (1987). For the SAHEX 
coatings characterized by separable amplitude transmission 
relations, the diagonal and circular polarization components, 
0 ( 2 ,*) and c(3,n>, are always zero. The normalized C vector for 
the separable amplitude transmission relations is 


t t l+d(0,2)i 2 + , d( l,0)+d( l,2)i + ... , 0, 0 ] 


where 


t . C (0, 0) and d(k,n) • c(k,n)/c(0,0) . 


. • r vector for coordinates other than the 

The Jones matrix and C vector 


s-p coordinates are 


obtained from the polarization rotation 

coordinates are rotated with 


operation. 


For example, the s-p 


A\ 


bv 0, the orientation of the plane 
respect to the x-y coordinates by 0 , 

-r mahrlx in x-*v coordinates JJ(x,y) i s 

of incidence. The Jones matrix in x y 

_ _ mat-rix in s-p coordinates JJ(sp) by the 

related to the Jones matrix in s P 

equation 

jj(x,y) - Rt-9) «<•'»> R(9) ’ 

The Taylor aeries coefficients for the Fresnel equations 
which govern an uncoated dielectric or metal surface have been 

determined for use in determining the instrumental polarisation 

of a conventional Cassegrain telescope. The notation a(s) 
a«p, will refer to either the reflected or transmitted amplitude^ 
transmission coefficient, while t(s), tip), rls) and rip) are 
used to refer unambiguously to the transmitted or reflected 
components. The Fresnel amplitude transmission equations are: 

t(syi) ■ ( 2 cos 1 sin i ' >/ (sin( i+1 ’ )) 

*’ = ( 2n cos i)/(n cos i + n'cos i’> 

t(p , U . ( 2 cos i sin i ' )/( sinl i+i ' ) cosli-iM) 

. ( 2n cos l)/l n’ cos i+ n cos i') 

r|s,i) - |-sin(i-i'))/( sinl i+i ' ) ) 

. , „ COS l - «' cos i')/< n cos i + n' cos i 1 ) 


r(p,i) - I tan I i-t ' ) )/( tan! 1+1 ' ) ) 

3 ( n' cos i - n cos i* )/( n* cos i +n cos i') • 


The Fresnel equations depend on the ratio of the indices, n and 
n', but not on the values of the refractive indices 
individually. This relative refractive index ratio is defined as 

N = n/n' 

The Fresnel equations are equally valid for real n, corresponding 
to transparent media, or complex n, corresponding to absorbing 
media or metals. 

The second order Taylor series expansions for the Fresnel 
amplitude coefficients about i=0 are: 

t(s,i) = ( 2 N)/(N+1) + i 2 (N(N-l) )/(N+l) , 

t(p, i ) = (2 N )/(N+l) + i 2 (N 2 (N-1))/(N+1) , 

r ( s , i ) * (N-l )/ (N+l ) - i 2 ( N-l )/ (N(N+1 ) ) , 

r(p,i) ( N-l ) / (N+l ) + i 2 (N-l )/ (N(N+1 ) ) • 

The direct method for calculating the Taylor series 
coefficients of a coating series given in the last section are 
impractical for multilayer coatings due to the complexity of 
calculating the partial derivatives of the appropriate amplitude 
transmission equations. The Taylor series coefficients can be 
obtained numerically from the s and p amplitude transmissions 
evaluated at a series of angles of incidence. An algorithm to 
sixth order has been given by Chipman (1987). The algorithm was 
used with the thin film design programs Filmstar and Films to 


obtain the Taylor series expansions of the transmitted and 
reflective amplitudes as a function of the angle of incidence for 
use in the polarization aberration calculations for SAMEX. 

CASCADED WEAK POLARIZERS 

In this section the Jones matrix describing the instrumental 
polarization for light propagating along a ray path through the 
SAMEX foreoptics is derived. Results are also given for the 
instrumental polarization associated with paraxial rays as 
functions of the Taylor series of the C vectors representing the 
optical interfaces. The notation used in this section is 
compiled in Table C-3 . 


Table C-3 . Notation for Section 



C 

c(k) 

d(k) 

i 

JJ 

k 

d 

q 

Q 

v 

r (k) 

a 

t 

p(k) 


d(k)°coeff relents rotated Into arbitrary plane of 

Normalized C vector components In s-p coordinates 

Artyle of incidence 

Jones matrix , , 

Pauli spin matrix index: 0,i,z,J 

Length of a ray segment 

Surface index 

Total number of surfaces 

Orientation of the plane of 

Absorption or polarization coefficient 

Pauli spin matrix 

Normal transmittance . 

Phase or retardance coefficient 


Subscript Ordering: k, 1, q' 


, _ i o i \ is the coefficient for: 
the vu" Polarization basis state, that is second 
orieO in the angle of incidence Taylor series, 
i 5 , for q=3 , the third interface. 



Consider an optical system with Q optical interfaces 
numbered in the order encountered from ,-l to 0. Ho symmetry 
regarding the optical configuration is assumed. Light propaga e 
along a specified ray path such as would be calculated by an 
optical ray trace calculation. At each interface 
polarization is introduced due to differences in the optica 
constants across the interface. In addition, polarization is 
associated with the ray path between interfaces due to optica 
active crystals, dlchroism, birefringence, gradient index 
materials or other polarizing mechanisms. But for the 
polarization analysis for SAHEX foreoptics, polarization 
associated with the optical path between interfaces -as zero 
Therefore only interface induced polarization is considered ere 


Homogeneous Optical Systems 

A homogeneous interface has optical properties independent 
of spatial coordinates on the interface, the dones matrices are 
functions only of the angle of incidence, plane of -idence, and 

optical properties of the interface media, JJ - JJ(l ' 
similarly C - C(i, 9 ,n,nM. the foreoptics and polarimeter 
sections of the SAMEX Kagnetograph are homogeneous optica 

systems. 

Likewise , a homogeneous medium has optical properties 

, . An anisotropic crystalline 

independent of spatial coordinates. 

e if it consists of a single crystal. The 
medium is homogeneous if it 

i *-h direction but not with position, 
refractive index varies with direct 


Radially Symmetric Systems o£ Lenses, Mirrors and Coatings 

The polarization properties of optical systems comprised of 
lenses, mirrors and "fine" coatings will be developed. A 
radially symmetric optical system has an axis of symmetry, the 
optical axis. It is assumed that the optical elements and 
materials used in transmission are highly transparent and 
nonpolarizing, as is usual in lenses. The polarization 
contribution from the path lengths through highly transparent 
elements is small relative to the polarization arising at the 
interfaces and is neglected. 

The polarization associated with ray paths near the optical_ 
axis, or in the paraxial regime, will be derived. For this 
paraxial development to be accurate, it is only necessa y 
the angles of incidence are small enough that the polarization 
associated with the interfaces is adequately approximated by a 
second order expansion of the C vector as a function of the angle 
of incidence. For an uncoated lens or mirror, this approximation 
is generally valid for i. < 30 degrees. Calculation of the fourth 
and higher order coefficients allows estimation of the accuracy 
of these second order equations. The paraxial region for this 
polarization analysis is typically orders of magnitude larger 
than the paraxial region of geometrical optics <the region where 
the fourth and higher order wavefront aberrations are 

negligible . ) 

Homogeneous and isotropic interfaces do not display 

-l • ■dcnoo There is only an amplitude and 
polarization at normal incidence. There 



Homogeneous and isotropic interfaces do not display 
polarization at normal incidence. There is only an amplitude and 
phase change which is represented by the complex number, t, the 
normal amplitude transmittance. An isotropic interface such as a 
lens, mirror or coating has a C vector Taylor series in s—p 
coordinates (9 * 0) of the form 

c(i ' 9 ) l 0 = o * t [1,0, 0,0] + i 2 t [d(0,2) ,d( 1,2) ,0,0] + 

+ i 4 t [d(0,4) ,d(l,4) ,0,0] + ... 

For an arbitrary orientation 9 of the plane of incidence, the C 
vector is i 

C(i,9) = t [1,0, 0,0] + i 2 t [c(0,2)t,c(l,2)t,c(2,2)t,0]- 

+ i 4 t [c(0,4)t,c(l,4)t,c(2,4)t,0] + ... 

where the c's are determined from the d's by a rotational change 
of basis. Since homogeneous and isotropic interfaces do not 
display circular retardance or circular polarization, a( 3) is not 
included to simplify the mathematics. 

The C Vector for a Paraxial Ray 

The SAMEX instrumental polarization will be analyzed by a 
paraxial optics development. Consider a paraxial ray path 
through an optical system from surfaces q=l to 0 with angles of 
incidence, i(q), and orientations of the plane of incidence, 

9 (q ) . The Jones vector associated with the axial ray (down the 




optical axis,) i«J>-° f° r a11 q ' U 

t(q-l) a( 0 ) . ■ • t C 2 ) a(0) t(l) c(0) =■ T <r (0) 


jj = t(q) <?(°) 


where 


T = 


Q 

it 


q = 1 


t(q) 


The complex amplitude transmittance down the axis. T, 

product of the normal incidence complex amplitude transmittance 

at each surface. (nter £ace q can be 

The Jones matrix associated -ith a ray a 

nf the interface Jones matrix 
expressed in terms of the expansion of the i ^ 

as 


. t(q) [e(0)+ i 2 , q Hc( 0 , 2 ,q)e( 0 ,tc(l, 2 .q)e(l)t=U, 2 .q)e( 2 ), 

MU ' V ' q) ‘ + 1 ^ (q )(c(0.4.q)e(0) + c(l»4.q)°(l^ + c(2»4.q)o(2)) + ...l 


, rix ass0 ciated with the entire paraxial ray 
The Jones matrix assoc 

path resulting from Keeping terms to second order at each^ ^ 
interface is (where X represents multiplication came 

next line) 


JJ = 


t(q) (o(0)+ i2(q)c(0,2,q)c(0) + c(l,2,q)ea)«(2-2'-^> 3(2,, * 

t(q .lHe(0 )+ i 2 ,q-U,c«0.2,-l,e<0Hc(1.2,-l,ea : d ;: ne,2, 


)x 


t(2)(a(0) + i 2 (2)(c(0,2,2)a(0) + c(l,2 / 2)a(l) + c(2,2,2)a(2))X 

t( 1 )(a (0) + i 2 (l)(c(0,2,l)a(0) + c(U2,l)a(l) + c(2,2,l)a(2)). 
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Associated with each interface are four terms. Carrying out all 
the multiplications leads to 4° terms, all in even powers in i. 
Collecting terms of equal power in i, there is one term at 
zero' th order and 3Q terms at second order. If i is assumed 
small, the large number of higher order terms are of diminishing 
importance. Collecting zero and second order terms in JJ Q and 
JJ 2 the expression for JJ is 

0 

JJq + JJ 2 * T j(0) + T <j( 0 ) E i 2 (q)c(0,2,q) 

q=i 

Q 2 

+ T a( 1 ) E i 2 (q)c( 1 , 2 ,q) 

q=i 
0 2 

+ T o(2) E i 2 (q)c(2,2,q) . 

q=l 

Since no polarization or retardance was assumed on axis, the 
contributions to the second order polarization for this ray are 
just sums of contributions from each surface. The multiplication 
taking place at second order for the a(l) term is of the form 

i 2 (q)a(l) a(0) a(0) a(0) . .. a(0) + a( 0 ) i 2 (q-l) a( 1 ) <y( 0 ) o( 0 ) . . . a( 0 ) + 

... + a( 0 ) o( 0) a( 0 ) a( 0 ) . . . i 2 ( q ) <j( 1 ) 

where the c-dependence is not shown explicity. 

This equation contains the useful result that, when no 
elements display polarization or retardance at normal incidence, 
as in the SAMEX foreoptics, there is no order dependence in the 



second order terms. Only one non-identity matrix term occurs in 
each second order matrix product. The second order polarisation 
associated with the paraxial ray path is obtained by a simple 
summation of second order polarization contributions at each 
intercept. Chipman (1987) gives a complete account of this 

derivation. 


Paraxial Optics Geometry 

The polarization aberrations for SAMEX are a description of 
the polarization behavior of an optical system expressed as an 
expansion about the center of the object and the center of the 
pupil. Thus it is appropriate and convenient to obtain the 
derivations from a paraxial ray trace; appropriate, because 
understanding the instrumental polarization near the center of 
the pupil and image is key to understanding instrumental 
polarization in general; convenient because the paraxial ray 

trace is linear, and thus easy to manipulate. 

The paraxial coordinate system used is a normalized nghU 

handed coordinate system. The z axis is the optical 

4 . . S vstem. Light initially travels 

rotationally symmetric optical system 

• „ „ Pinure C-l shows the notation, 
in the direction of increasing z. Figure c 

For a rotationally symmetric system, the object can be 

„ „ h . „ „ is without loss of generality. The obje ct, 

located on the y axis wituwu*. 

, w n,af h = 0 in the center of the 
coordinate H is normalized such tha 

• \ anH u s l at the nominal edge of the 

field (on the optical axis) and H 


field o£ view* 

„ <? t-rik.es the entrance pupil is 

The location where a ray strides 



r l Paraxial coordinate system. The paraxial system is 

a normalized right-handed «lfo^lca5%m^ -light 

SUE ^“-av Is n tb, a dir. ^lon o< ^-eaaln, z ^ « 

an optical system are =* arae * e £*®?. * * ' Uze d object coordinate, 
object and entrance pupil. H Is the normalize J__ e ln the 

p is the normalized pupil radius, and ♦ z p The normalized 

pupH meas'QTfed counterclockwise from the y axis. The norman 

Cartesian coordinates in the pupil are x and y. 

marginal rays are also shown. 
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specified by the polar pupil coordinates p and ±. p is 
normalized such that at the edge of a circular pupil p » 

The angle * is defined here as it is in much of geometric optics, 
and in defiance to most analytical geometry, as being zero on the 
•y axis' and increasing counterclockwise. Normalized Cartesian 
pupil coordinates x and y can be used. They are defined as: 
x 3 p sin( $) and y - p cos($). 

Expressions for the angle of incidence i and the orientation 
of the plane of incidence ♦ of a ray at a given surface q will be 
expressed in terms of the marginal i(m,q) and .chief ray i(c,q) 
angles of incidence at that surface. Details of the derivation 
may be found in Chipman ( 1987, section D) . Note, however, for a 
radially symmetric system the angle of incidence should be a - 
function of H 2 , p 2 , and Hpcos* since the function should be 
invariant to rotation of the system about the optical axis and 
must give the same result when -x is substituted for +x where 

x - pcos*. 

Assume that a paraxial ray trace has been performed for a- 
specific system and that i(m,q) and i(c,q) have been 
calculated. A ray from normalized object coordinate H which 
passes through pupil coordinates p and * will have an angle of 
incidence i(q> and orientation of the plane of incidence 9(q) at 
surface q equal to: 


i (q ) 


= [H 2 i 2 ( c ,q ) 


+ 2Hp cos(-fr) i(c,q) i(m,q) + p 2 i 2 (m,q)] ^ 


= sin" 1 [ p sin( <fr) i(m)/| i 1 1 
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Figure C-2 shows the paraxial angle and plane o£ incidence 
for three field angles. The magnitude of the angle of incidence 
is represented by the length of the lines. The orientation of 
the lines corresponds to the orientation of the plane of 
incidence. Note that off axis, the pattern is a shifted version 
of the on axis pattern. The incidence angle is given by the 
pupil coordinate (p,«) and the image coordinate (H> since these 
coodinate define the optical path of a single ray. 


POLARIZATION ABERRATIONS 
Introduction 

This section derives the polarization aberrations for SAMEX - 
as a Taylor series description of the instrumental polarization 
associated with paraxial rays through the foreoptics. 

Polarization aberration is a description of the polarization 
behavior of an optical system expressed by a expansion about the 
center of 'the object and center of the pupil. Table C-4 gives an 
overview of polarization aberration theory. The results are 
obtained in a form very similar to the wavefront aberrations. In 
particular, terms closely related to defocus, tilt, piston error 
as well as the Seidel and higher order aberrations can be 
associated with all eight of the basis Jones matrices. Since 
polarization effects are typically orders of magnitude smaller 
than wavefront effects, fewer terms are needed for a sufficient 
description. A method of calculating aberration coefficrents for 

is developed in the next section. 


specific systems 
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Polarization Aberration Theory Outline 



(1) Jones Vector 

The Jones vector is a complex 
2 — component vector describing 
the electromagnetic field. 


(2) Paraxial Ray Trace 

The angle of incidence, i, of any 
ray can be written in terms of the , 
angle of incidence of the chief 
and marginal rays, i(c) and i(m). 

(3) Coating Calculation 

Thin film design programs 
calculate the amplitude 
transmission coefficients for s 
and p light, t(s,i) and t(p,i). 

(4) Taylor Series Representation 

The amplitude transmission 
coefficients, t(s,i) and t(p,i)» 
are transformed to with a Taylor 
series expansion about normal 
inc idence . 

(5) Jones Matrix for an Interface 

The second order Jones matrix for 
an interface can be written in terms 
of three Polarization Aberration 
terms represented by coefficients. 

(6) Aberrations Sum for System 

The second order polarization 
aberrations for the system is 
the sum of the aberration 
contributions of each interface. 


i ( c ) , i(m) 


t(s,i)r t(p,i) 


t(s,i) * t Q (l + tCs,2)i 2 + ... 
t(p,U = t Q (l + t(p,2)i 2 + ... 


P(l,0,2,2,q) , P( 1,1,1, l»q) * 

P( 1 , 2 , 0 , 0 ,q ) 


P( 1 , 0 , 2 , 2 ) = q 

P(l, 1,1,1) =* g 
P(l,2,0,0) = q 


P(l,0,2,2,q) 
P( 1 , 1 , 1 , 1 ,q) 
P( l,2,0,0,q) 


(7) Polarization Accuracy 

Summarize the performance of the system 
with a single number, the polarization 
accuracy Up). This is the maximum 
fraction of light coupled into the . 
orthogonal polarization state* ^ ^ , 
occurs at the edge of the field of view. 




P 


3 


The Polarization Aberration Expansion 

The wavefront polynomial expansion describes the variation of 
the optical path difference through an optical system as a function 
of ray coordinates. A closely related expansion will be presented 
for all four basis polarization matrices «(0>, .ID. v<2), °< 3) - 
The polarization aberration expansion for radially symmetric 
systems uses a very similar polynomial expansion to describe all 
eight basis polarization vectors. The principal difference is a 
modified form lor the linear polarization and linear retardance 
terms since these involve both a magnitude and an orientat 
The eight forms of polarization behavior can be 
characterized by four complex numbers, for example, the four 
elements of either the Jones matrix or the C vector. Me 
introduce a new set of complex parameters, the complex 
polarization aberration coefficients which gives a description of 
the polarization behavior of an optical interface. It should be 
emphasized that the amplitude and phase of the coefficients are 
generally unrelated. They refer to different aspects of the 
instrumental polarization. The amplitude part of the coefficient 
describes amplitude and polarization effects while the phase part 
describes phase and retardance. When necessary the amplitude and 
phase of the polarization aberration coefficient, P, will be 
denoted by A and $ where, 


P a A exp( t $) 


#•» 



The complex aberration coefficient is written to contain 
polarisation, amplitude, effects, retardance, and phase effects 
in a single term. The following polarisation aberration 
expansion for seroth and second order is used. The subscripts 
are defined as follows : 

P(k,u,v,w) = A(k,u,v,w) exp(T*(k,u,v,w) ) , 

where: k is the type of polarization behavior, 

u is the order of the H dependence, H , 

v is the order of the p dependence, P V , and, 

w is the order of the <> dependence, (cos(*)) w . 

The indices u, v, and w are used exactly as they are for the 

wavefront aberrations, as shown in the next section. 

The polarization aberration expansion of the Jones matrix 

for the SAMEX foreoptics is 

3 ' ‘ 

JJ(H, p ,4>) * £ c(k,H, ,r,p) <j(k) 

k=0 

= \ Z E Z P(k,u,v,w) H u p v cos w ($) cr ( k. ) . 

k=0 u v w 

Here we define the aberration coefficients as an expansion of the 
Jones matrix in terms of the ray coordinates p, <fr, and H. 

The C vector coefficients in this polarization aberration 

expansion are: 


3 . 


Amplitude and Phase term: 


c ( 0 / H t p t ♦ ) - 


A(0, 0,0,0) + A(0,2,0,0)H 2 
+ A(0, l,lrD H p cosU) + A( 0 , 0,2,0) 

+ t( *(0,0, 0,0) + ♦ ( 0 , 2 , 0 , 0 ) H 2 2 

+ *(0,1,1,1)H P cosU) + MO, 0,2,0) P 


Linear 


Polarization and Retardance terms: 


c ( 1 1 H t P f ♦ ) 


A(l, 0,0,0) + A( 1,2, 0,0) H 

* H p (Ml.lrl.U COSU ’ ' M2,1 ' 1 ' 1) 

* p 2 (M1, 0,2,0) oos(2*) - M2,0,2,0) 

+ T (*d, 0,0,0) + *(1,2, 0,0) H 

H p (*(1,1, l.D cos( *) ' * 12 ' 1 ' 1 ' 11 

p 2 ( *(l, 0 , 2 , 0 ) cost 2*) - »(2, 0 , 2 , 0 ) 


Diagonal" Polarization and Retardance terms: 


c(2fHf Pr 


A(2, 0,0,0) + A(2,2,0,0) H 2 + 

+ H p (A(2, 1,1,1) cos(*) + A(l, 1,1, 
+ p 2 (A(2, 0,2,0) cos ( 2 *) + A(l, 0,2, 
+ T (*(2,0, 0,0) + *(2,2, 0,0) H 2 + 

H p (*(2,1, 1,D cosU) + Hl,l,^ 
p 2 ( *( 2 , 0 , 2,0) cos ( 2 *) - Ml, 0,2. 


Circular 


polarization and Retardance 


terms : 


) , 


sin ( * ) ) 
sin ( 2 * ) ) 

sin( *) ) + 
sin( 2 *) ) ) • 


1) sinU)) 

0) sin (2*) ) 

1) sin( *) ) + 

,0) sin(2*) )) 



c(3,H,p,*) 


= A(3, 0,0,0) + A(3,2,0,0)H 2 + 

2 

+ A(3, 1,1,1) H P cos(<J») + A(3, 0,2,0) p 
+ T (*(3,0, 0,0) + *(3,2,0,0)H 2 

2 

+ * (3,1,1,1)H p cos(*) + ♦( 3 , 0 , 2 , 0 ) p )/ 

There are thirty two terms in this polarization aberrat' 
expansion to second order arising from four terms in each of the 
eight degrees of freedom of the Jones matrix. The terms y 
grouped as follows: 


A( 0 ,u , v ,w) 

Amplitude terms 

A( l,u,v,w) 

Linear polarization terms 

A( 2,u,v,w) 

Diagonal polarization terms 

A(3,u,v,w) 

Circular polarization terms 

*(0,u,v,w) 

Wavefront or phase terms 

*(l,u,v,w) 

Linear retardance terms 

* ( 2,u, v,w) 

Diagonal retardance terms 

* ( 3 ,u , v,w) 

Circular retardance terms 

P(k, 0,0,0) 

"Constant Piston" terms 

P(k, 2,0,0) 

"Quadratic Piston" terms 

P(k, 1,1,1) 

"Tilt" terms 

P(k, 0,2,0) 

"Scalar Defocus" terms 

p (k , 0 , 2 , 2 ) 

"Vector Defocus" terms 


The names of the wavefront aberrations: piston, quadratic 

used here in an extended sense, to 


piston, defocus and tilt, are 



describe variations of components of the Jones vector which share 
the same functional dependences as the wavefront aberrations. 
Defocus is a p 2 variation of a parameter. Tilt is H p cosU) 
variation. Quadratic piston is H 2 variation. Thus, "amplitude 
defocus" means a p 2 amplitude variation. Likewise the "circular 
retardance tilt” is a the H p cos(*) circular retardance 
variation, and so on. 

This polarization aberration expansion is an equation which 
describes all possible second order variations of the Jones 
matrix, just as the second order wavefront aberration expansion 
spans the set of all second order wavefront variations. Thus the 
polarization aberration expansion characterizes quadratic 
variations of all forms of wavefront, amplitude, polarization and 

retardance . 

This polarization aberration expansion is a summation of 
terms in the different Pauli spin matrix components, not a 
product. Thus the four C vector elements can be pictured as 
acting in parallel, almost side by side in the aperture, but not 
in series. Each term describes an amount of a particular form of 
polarization, independent of the other contributions. 

An "aberration term" is to be considered as containing all 
the algebraic terms in the expansion with the same coefficient. 
Most of the coefficients occur only once and the aberration term 
contains only one algebraic term. The exceptions are the terms, 
A(l, 1,1,1), *(1,1, 1,1)/ A( 1 , 0 , 2 , 2 ) , *(1,0, 2, 2), A(2, 1,1,1), 

*(2,1, 1,1), A( 2 , 0 , 2 , 2 ) , and *(2,0, 2, 2). These aberration terms 
have components both along the axes and at 45 degrees. 


With SAMEX, the principal concerns are with the linear 
piston, linear tilt and linear deEocus terms, both in 
polarization and retardance. These are going to be the largest 
terms present which corrupt the incident polarization state. 

These values are given in the Table 19 for SAMEX. 

For a detail discussion of the physical meaning of the 
polarization aberration coefficients see Chipman (1987), however 
a discussion of the orgin of tilt and piston, P(0,2,0), P(l,l#l)» 
and P ( 2 , 0 , 0 ) terms is included below. 

A distinction is made between scalar and vector 
aberrations. The wavefront aberrations are scalar aberrations, 
single valued functions of object and pupil coordinates. The 
linear polarization and linear retardance aberrations are vector 
aberrations since a magnitude and orientation is associated with 
these at each point. Amplitude, circular polarization and 
circular retardance aberrations are scalar since they are single 

valued and range positive and negative. 

Figure C-3 (top) shows the chief and limiting rays at an 

interface for objects on axis and at the edge of the field of 
view. Figure C-3 (bottom) is a plot of the value of the angle of 
incidence along the y axis as a function of p. Tilt terms 
naturally occur because as the object point moves off axis, the 
angle of incidence increases at one edge of the beam and 
decreases at the other edge. Tilt contains the first order 
portion of . this correction. 

Figure C-4 shows the off-axis angle of incidence squared and 
the decomposition of this into defocus, tilt and piston terms. 


OFF-AXIS 




ON-AXIS 


ANGLE OF 
INCIDENCE 



COORDINATE 


Figtire Q Angles of incidence for objects on and off axis. I 
e top frame, the chief and limiting rays at an interface are 
shown for objects on the optical axis and at the edge f the 
field of view, in the lower part of the figure, the angle of 

f °f lncident alon 9 the y axis (in the pa?axial 

system) is plotted as a function of the normalized radius o 
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Off-axis 
Angle of 
Incidence 
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Defocus 


2ax 

Tilt 


2 

a 


Piston 


igure f tlllT angle°of i^cidenle is 

S SSofitK ^“diU tilt and piston 


of incidence. 

shown along with 
terms . 


These terms are required 
vertex is located at an 


to describe a quadratic var 
arbitrary position on the y 


iation whose 
axis because 


x 2 ♦ 


(y-a) 


2 = x 2 + y 2 _ 2y a + a' 


this case since a is a linear function of «, 


= k H 


then the 


quadratic polarisation variation becoa.es 


c(x 2 + (y - kH > 2 ) - c(x 2 *y 2 ) - 2C y k H ♦ C K 2 H 2 

* p( 0 , 2,0 ) 0 2 +P(1,1.1)2H 0 cosU) * P( 2,0,0) 


where P(0,2,0), PfU.U, and P(2,0,0> are the defocus-, tilt and 

quadratic piston aberration coefficients and the 
index is not shown. Tilt and piston tenns arise nature y 
decentered defocus. Similarly, the fourth order -efront 
aberrations co.a, astiq.atis., field curvature and dis or 
ari se naturally fro. decentered spherical aberratron, a • 



CALCULATION OF ABERRATION COEFFICIENTS 

The method used to calculate the second order polarization 
aberration coefficients for the SAME* foreoptics given the C 
vector power series for each interface is detailed in this 

section . 

Single Surface Aberrations for Amplitude Transmittance Relations 
For homogeneous and isotropic interfaces characterized by 
amplitude transmittance relations, such as lenses, mirrors and 
thin film coatings, the polarization aberrations at a inte 
simplify considerably. At these interfaces the Fresnel equations 
and related thin film equations are separable into s and p 
components, so the Jones matrices representing the interface in 
a-p coordinates are diagonal. The off-diagonal terms, diagonal 
polarization «(2) and circular polarization ,«(3) , are not 
present. Further, with isotropic media, the s and p amplitude 
transmission coefficients at normal incidence must be equal. ' 
Thus the amplitude transmission functions can be expanded as: 


a ( s 1 1 ) 


a(0) ( 1 ♦ a(s,2)i 2 exp(T(d(0) ♦ d(s,2)i 

t (l+(a(s,2) * T d(S,2))i 2 +...) . 


4 ^ + . • • ) 


a(p,i> = a < 0 ) ( 1 + a(p,2)i 2 ♦. . . ) exp(T(d( 0) + d( 2,p) 

* t (1 + a(p,2) + T d(p,2)i +••♦) 


t2 +. . . ) 


where : 



t 


a( 0) exp ( T d( 0 ) ) . 


The s-p coordinate Jones matrix expansion to second order is 


JJ( i) 


a ( s , i ) 0 

0 a(p, i ) 


= t 


( a ( 0 ) ( 1 + A (2) 


7D(2))i 2 ) + a(l)(a(2) + id(2))i ) 


where 


^(2) = ( a ( s , 2 ) +a (p, 2 ) )/2 , 


a(2 ) » (a(s,2)-a(p,2))/2. 


D(2) 


= (d(s, 2 )+d(p, 2 ) )/2, 


d(2 ) = (d(s,2)-d(p,2))/2. 

The s-p coordinate C vector Taylor series expansion to second 


order is 


C(i) - t (1/0, 0,0) 


+ i 2 t (A(2)+TD(2), a(2)+Td(2), 0, 0) 


The x-y coordinate C vector Taylor series for orientation of the 
plane of incidence 0 is 



0(1,0) » t(l, 0,0,0) + i 2 t (A( 2 ) +TD( 2 ) , 

(a ( 2 )+ id ( 2 ) ) cos (29), (a(2) + id( 2 ))sin( 29 ) 

The normal -incidence polarization aberration terms (the constant 
piston terms) are zero: 


P(l, 0,0,0) 


P(2, 0,0,0). = P(3, 0,0,0) = 0 . 


There 


is no polarization or retardance on axis 


only the 


amplitude and phase transmission factor t. 

Ail terms for off-axis diagonal and circular polarization 
are zero: 


p (2,u,v,w) = P( 3 ,u , v,w ) = 0 . 

Thus, the single surface C vector in paraxial coordinates is 
obtained by substituting the paraxial representation of i(H,„,*) 
and v(H,p,$) into C(i,v) yielding 

c(0,H,p,^) » t + t d(0,2){ H 2 i 2 (c) 

+ 2H p cos(*) i(c)i(m) + p 2 i 2 ( m )] 

c(l,H,p,*) = t d ( 1 , 2 ) ( h 2 i 2 ( C ) 

+ 2H p cos ( $ ) i(c)i(m) + p 2 cos(2$) i 2 (m)] 

c(2,H,p,^) = t d ( 1 , 2 ) [ 2H p sin(<fr) i(c)i(m) 

+ P 2 sin ( 2 $) i 2 ( m ) ] 


H2 


c(3,H,p,$) » 0. 

Since there Is no diagonal polarization, the only contributions 
to c<2) arises from the rotation of linear polarization from the 
S-P coordinates into the x-y coordinates. 

Polarization Aberration Coefficients for Systems 

Since the polarization aberrations are only being evalua 
to second order in -the angle of incidence, the difference between 
spheres, parabolas, conics or other radially symmetric aspherics 
does not occur at this order. The relevant shape parameter here 
is only the vertex radius of curvature. The angle and plane of 
incidence differences between these types of interfaces are the 
same at second order but will differ at fourth order and higher. 

For surfaces q-1 to 0, each surface is characterized by 
three complex parameters from the normalized C vector expansion: 


d (0,0) = t(q), d(0,2,q), and d(l,2,q) 



The single surface 


polarization aberration coefficients are 


P(0,0,0,0,q) 

P(0,2,0,0,q) 

P(o,i,i»i#q) 

P(0,0,2,0,q) 
P( l,2,0,0,q) 
P( 1,1,1, l,q) 
P(l,0,2,2,q) 


* t (q ) 

= t ( q ) d ( 0 , 2 ,q ) i 2 (c) 

= 2t(q) d ( 0 , 2 ,q) i(c)i(m) 
= t(q) d( 0 , 2 ,q) i 2 (m) 

= t (q ) d ( 1 , 2 ,q) .i 2 (c) 

= 2t(q) d( 1 , 2 ,q ) i(c)i(m) 
= t (q) d ( 1 , 2 ,q) i 2 (m). 


The polarization aberration coefficients for the system are 
calculated by chain multiplying the single surface polarization 
aberration expressions and keeping terms to second order in H and 
p. The zero and second order Jones matrices for the g'th 

interface are: 


jJ 0 (q,H,»,*> - P(0,0,0,0,q) a(0) - d(O,0,q) «(0) = t(q) olO) . 

JJ a (q.H.».*> - .(0) (H* P(0.2.0.0.q» * 2H p cost 

+ p 2 P ( 0 , 0 , 2 , 0 ,q) 1 


+ o( l) [ H 2 P(l,2,0,0,q) + 2H p cos U ) P ( 1 , 1 , 1 , 1 rQ > 

+ p 2 cos (2^)P(l#0»2,2,q) 1 

+ p( 2 ) t 2H p sin(*)Pd,l.l,l.q) * p 2 sin(2*)P(l,0,2,2,q> 1 



M. 


= a ( 0 ) t(q)d(0, 2) ( H 2 i 2 (c) + 2H p cos ( * ) i (c ) i (m ) 

+ p 2 i (m) 2 .] 

+ a ( 1 ) t(q)d( 1,2) [H 2 i(c) 2 + 2H p cos ( <►) i ( c ) i (m ) + p 2 cos(2$) i 2 ( 

1 

-t- a (2) t ( q ) d ( 1 f 2 ) [ 2H p sin( ♦) i(c) i (m) + p 2 sin( 2 *) i 2 (m) ] 

Multiplication of the single surface Jones matrices yields 

JJ(H,r,p) = ' JJ ( q , H , p , $ ) “ 

q=Q,-l 

= “ ! [ JJg ( q ) + JJ 2 (q,H,p,$) ]. 

Since JJg(q) is a constant function, independent of H, p, and 
the ( H , p , $ ) dependence can be dropped. This expression contains 
2^ terms. The order of a term is the sum of the powers of H 
and p, H u p v cos($) w > i.e. order = u + v 

There is one first order term and Q second order terms. 


The zero order Jones matrix is 



Ml 


JJ 0 - q"o,-l W 0 ( 5' B "'* n = q-1 t(q) = T ' 

the system amplitude transmittance. The second order Jones 
matrix is greatly simplified since, for isotropic surfaces, all 
zeroth order Jones matrices are a constant times the identity 
matrix o(0). The second order only contains products which 
contain a single second order term.’ The second order Jones 

matrix is 

JJ 2 (H,p, 0) = q=l,Q * 

At second order the weakly polarizing isotropic interfaces do -not 
display order dependence. The product of any two second order 
terns is fourth order. The order dependence enters at fourth and 
higher order. Second order is a simple sum of polarization 
contributions. Collecting the piston, tilt and defocus terms 
from the second order Jones matrix provides the coefficients for 
the system polarization aberration expansion to second order: 


P(0, 0,0,0) 
P( 0,2, 0,0) 
P(0, 1,1,1) 
P(0, 0,2,0) 
P ( 1 , 2 , 0 , 0 ) 

P(l, 1,1,1) 
P( 1,0, 2, 2) 


T 

T 

2T 

T 

T 

2T 

T 


, 


E 


q=l ,Q 
E 


q-l,Q 

E 


q=l ,Q 
E 

q*l,Q 


E 


q=l,Q 

E 

q=l,Q 


d ( 0 , 2 ,q) 
d ( 0 , 2 ,q) 
d ( 0 , 2 ,q) 
d( 1 , 2,q) 
d ( 1 , 2 ,q) 
d{ 1 , 2 ,q) 


i 2 (c) , 
i(c)i(m) 
i 2 (m) , 
i 2 (c) , 
i (c) i (m) 
i 2 (m) 
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The other three zero order =oe££icients and the other six second 
order coefficients are all zero: 


P( 1,0, 0,0) 
P(2,2,0,0) 

P ( 3 , 2 , 0 , 0 ) 


= P(2, 0,0,0) = P(3, 0,0,0) = 0 

» P(2, 1,1,1) = P(2, 0,2,0) = 0 

S P(3, 1,1,1) = P(3, 0,2,0) - 0 


The amplitude and polarization coefficients are the real 
parts of the P coefficients 


A(k 


,u,v,w) = Re(P(k,u,v,w)) . 


The retardation 


coefficients are the imaginary parts 


*(k,u,v,wT 


Im(P(k,u,v,w) ) . 


The polarization aberration coefficients are calculated for the 
foreoptics from the paraxial geometry and the normalized C- 
vectors for the coatings and interfaces. 


Polarization Accuracy 

The Jones vector gives the amplitude of the electric field 
and the square of the amplitude gives the intensity of the 
components. With the Jones matrix one is able to calculate the 
polarization effects for the optical system. From above the 
second order the Jones matrix which gives the linear polarization 

and linear retardence is 


JJ 2 (H,P,$) 3 a( 1 ) [H 2 P ( 1 , 2 , 0 ,0 ) + 2 H 0 cost*) 911 , 1 , 1 , 1 ) 

+ p 2 cos ( 2 $) P(l,0,2,2) ] . 

+ a( 2 ) [ ( 2Hp s in ( $ ) P(l, 1,1,1) + p 2 sin( 2 4 . ) P ( 1 ,0,2,2) 1 

The on axis linear polarization and linear retardance of the 
SAMEX foroptics are characterized by the term linear defocus. 

The instrumental polarization function JJ(H,p,$) for linear 
defocus is 

JJ(H,p,$) = t[<Jq + P( 1 ,0,2,2) p 2 (o^cos2$ - 02 s ^ n 2 'M ~~ 

- tU 0 ♦ ( a 10 22 + j«1022 )|,2l ‘'l CO32 * * J 2 sin2 * )1 ’ 

Here, t is the amplitude transmittance of the system down the 

optical axis. It ( x ) describes the polarization independent 

reflection and absorption losses associated with the ray down the 

optical axis at normal incidence at all interfaces. P(l,0,2,2) 

describes the linear polarization (a 1022 ) and linear retardance 

(5 ) associated with the marginal ray. 

1 u ^ ^ 

Maximum coupling occurs when the incident light is 
circularly polarized, since circularly polarized light can always 
be decomposed into equal components of J r and everywhere in 
the pupil. The coupling is zero in the center of the pupil 
(where the polarization and retardance vanishes) and increases to 

a maximum coupling fraction of 


n 


L C/itiax 


( H r 1 r ♦) “ |P(i'0,2,2)| - a 1022 * ^1022 


>t the edge of the pupU. The net fraction of incident circular 
polarized light coupled into the orthogonal circular polarized 
3ta te is given by the integral over the pupil of 


I c - 1 1 1 2 J 2 d* / pdp|p 2 rP(i ,°' 2 » 2 M 


0 o 
2 


_ tP ( 1 , Q , 2 , 2 ) 

5 - 

For incident linear or illiptical polarized light, the fraction 
of coupled intensity is less because the light is not compose 

equal fractions of eigenstates. 

The coupling is minimum for incident linear polarize 

ehich will be in one of the elgenpolarlzations along one axis in 
the pupil and in the orthogonal eigenpolarization along the 
orthogonal axis. The fraction of coupled energy will be 
calculated assuming an incdient polarization state of horizonta 
linear polarized light H for calculational simplicity, 
fraction is coupled for any linear polarized incident state, 
polarization state transmitted by an optical system described by 
linear polarization defocus for H in is 


J(H,Pr ♦> 


„ t( H ♦ P( 1 , 0 ,2,2) P 2 (S c° a 2 * - V stn U)) ' 


The fraction of incident H light couplied into V light is equal 



tP( 1,0,2,2) 


This is the mininmum fraction of energy coupled by linear defocus 
aberrations. Since any elliptically polarized incident beam can- 
be written as a sum of linear and circularly polarized light, the 
coupling fraction for arbitrarily polarized light lies in the 

range. . 

| tP( 1,0,2,2) 1 2 < j < |TP(l,0,2,2J_[ 2 t 
6 c 

When unpolarized of circularly polarized light is incident 
on the optical system, maximum polarization coupling occurs. 

This maximum is the polarization accurracy of a system and is 
calculated by the coupling integral^ I c , m ax' 


I - - / / 

c,max it 1 


pupil I 


c 2 (l) + c 2 ( 2) + c 2 ( 3) | pdpd$ . 


This integral must be evaluated numerically except for some 
special cases (see section 3.5). However, an analytic upperbound 
on the polarization accurracy can be easily established by using 
the triangle inequality 


c,max 


t S ! 


pupil 


[( (Re[c(l)] ) 2 + (Re [c( 2)] 2 


( Re ( c ( 3 ] ) 2 ) 2 


+( (Im(c( 1)] ) 2 + 


( Im [c ( 2 ) ] ) 2 + ( Im [c ( 3 ) ] ) ) 1 pdpd<t» 


We can evaluate the upper bound on the polarization accurracy 
using the second order polarization aberration expansion 
coefficients. The following integration is for the real part of 
the aberration coefficients. a iuvw 


. Uli / 2 * d, / l pd„[(a 1200 H 2 . a uu H» cos, + a 1022 » 2 cos2») 


2 2 

+ (a nll Hpsin^ + a 1Q22 P sin2*) ] 


2 h 2 a 2 

? 7 4 a llll H . 1022 

ti 2 u 2 200 h 4 + — i— + — r 


i . 


The corresponding integral 
substituting the imaginary 
coefficients, 5 i UVW/ ^ or a 


t. • is obtained by 

' 1 imaginary' 

part of the polarization aberration 
in the above expression. The 

luvw 


upperbound on the polarzation coupling is 


I Crmax * real 


< + ^imaginary 


= | T | 2 t a 1200 


„ . «? in H 2 6 io22 

h 4 + !miL ♦ non + 5 2 2oqH 4 + ♦ 


H 2 a 2 


,w |P(l.l,l,l)lV A lP(l>Q>2>2)_L 1t 

= |t| 2 [|P(1,2,0,0)| Z H 2 ~ 

This is the second order upper bound' to the polarisation coupling 

for systems of weak polarizers. 

Hence the average effect over the image and pupil can be 

obtain by integrating over a, and *. The polarisation 
accuracy, V defined as the ma*^ faction of light 
(intensity) which can be coupled into orthogonal polarisation 
state. The incident polarised state is given by the. Jones 
vector, 3. This vector is effectively rotated by the optical 
system and the rotation is given by the Jones matrix, JJ. The 
amount of- polarisation along the orthogonal state of 
polarisation, »•'. of the incident polarisation state is given by 

the projection of JJ(^) into J t 


JJ( 3) *3' . 


This value is given by the square of the second order Jones 
matrix and is given in terms of the polarisation aberration 
coefficients. (For the incident light in the polarisation 
state I p, then the result of instrumental polarisation is to 
couple the polarisation into the orthogonal state The amount 


of coupling is given by 


I P 
o r 


r o (i -y p + Vo q 


where A p » polarization 
polarization accuracy is 
polarization abberation 


accuracy.) From the above results, the 
given in terms of the second order 
coefficients P(l,0,2,2), P(,1»1»1»U» and 


P(l,2,0,0) : 


n 

a = - / I JJ-(H = l,p,$) | pdpd^ 

P TT J 1 * 


or 


2 (1,2, 0,0) + ip 2 (1,1, 1,1) + ^P 2 (1,0,2,21 


4p . 2P 


where the integrations have been carried out tor the squared 
terms and estimated tor the crossed terms. For the SA«EX 
magnetograph design given herein the polarization accuracy value 

of 


A p < 1.4 x 10 


u obtained tor specially selected optical coatings was 
obtained. These second order aberration coefficients are giv 
in section II-4. The second order coefficients are sufficients 
since the next order that contributes is the forth «der. The 
polarization effects would be on the order of (a p ) 



